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ABSTRACT 

We present a mechanism for chondrules to stick together by means of compaction of a porous dust 
rim they sweep up as they move through the dusty nebula gas. It is shown that dust aggregates 
formed out of micron-sized grains stick to chondrules, forming a porous dust rim. When chondrules 
collide, this dust can be compacted by means of rolling motions within the porous dust layer. This 
mechanism dissipates the collisional energy, compacting the rim and allowing chondrules to stick. The 
structure of the obtained chondrule-dust agglomerates (referred to as compounds) then consists of 
three phases: chondrules, porous dust, and dust that has been compacted by collisions. Subsequently, 
these compounds accrete their own dust and collide with other compounds. The evolution of the 
compound size distribution and the relative importance of the phases is calculated by a Monte Carlo 
code. Growth ends, and a simulation is terminated when all the dust in the compounds has been 
compacted. Numerous runs are performed, reflecting the uncertainty in the physical conditions at 
the chondrule formation time. It is found that compounds can grow by 1-2 orders of magnitudes 
in radius, upto dm-sizes when turbulence levels are low. However, relative velocities associated with 
radial drift form a barrier for further growth. Earlier findings that the dust sweep-up by chondrules is 
proportional to their sizes are confirmed. We contrast two scenarios regarding how this dust evolved 
further towards the densely packed rims seen in chondrites. 

Subject headings: solar system: formation — planetary systems: formation and protoplanetary disk 



1. INTRODUCTION 

Protoplanetary nebulae have been studied i n increasing 
detail from visual to mi crowave wavelengths (Mey er et alj 
120071: [Watson et al. 2007), and hundreds of extrasolar plan- 
etary systems have been discovered, but the 'primary ac- 
cretion' stage of the planetary formation process - that 
which leads from interstellar grains to planetesimals large 
enough to decouple from the nebula gas (asteroid- and 
comet-nucleus size objects) - remains obscure. In this par- 
ticle size range, coupled particle-gas dynamics dominate s 
the evolution, as review e d rece ntly by ICuzzi et alj d2005[) . 
ICuzzi & Weidenschillingl (120061) . and lDominik et all 00071) . 
The main processes which have been hypothesized for pri- 
mary accretion include i) incremental growth by sticking of 
small grains to each other and to larger particles, ii) various 
kinds of instabilities occurring in a particle-rich midplane re- 
gion, and Hi) formation of planetesimals from dense zones of 
particles that form in turbulence due to vorticity or pressure 
effects. 

A critical but unknown nebula property in this stage is 
whether turbulence is present, and if so, wh at its intensity is 
(IStone et al.ll2000HGammie & Johnsonl2005l) . If the nebula is 
nonturbulent, particles of all sizes can settle into a dense layer 
near the midplane where incremental growth is fairly robust 
for expected, but still poorly known and therefore somewhat 
ad-hoc, sticking properties (Weidenschillin g & Cuzzilfl993l : 



ICuzzi et alJfl99l IWeidenschillindfT997l l2000l 12004). This 
is because the dense particle layer drives the entrained gas 
to corotate, and relative velocities between equal-size parti- 
cles would largely vanish. The high local mass density en- 
sures that growth is rapid (|Weidenschillind l2()O0f) - perhaps 
too rapid (ICuzzi et alJ 120051 ). Various instabilities in such 
a layer, mostly gravitational, have been studied for decades 
( Gold reich & Wardll973t lSekivalfT998l: lYoudin & Shull2002T: 
lYoudin & Goodmanl2005l) . but these are precluded if the neb - 
ula is even weakly turbulent (ICuzzi & W eidenschilling 20061) . 

Astronomical and planetary observations seem to 
be most naturally reconc iled wit h turbulent nebu- 
lae (iDullemond & Dominikl l2005b ICuzzi et al.1 I2QQ1 



Electronic address: ormel@astro.rug.nl 

Electronic address: jcuzzi (iSmail.arc.nasa.gov, atielens@mail.arc.nasa.gov 



Brownlee etail l2006t IZolensky et al] l2006t ICiesla & Cuzzil 
20071) . A number of studies indicate that turbulence excites 
meter-size particles to relati ve velocities at which they 
probably disrupt each other dWei denschilling 1988; Benz 
2000; Sir onoll2004tlLangkowski et alj|2007l) . posing a barrier 
to further growth. However, some recent studies suggest that 
turbulence itself can concentrate particles of different sizes, 
in different ways, and trig ger rapid p lanetesimal formation 
(CuzzietaU 1200 ll 120071; [Rice et alj 120061 Uohansen et all 
120071) . Thus, in spite of the ongoing uncertainty i n just 
how turbulence may be maintained (IStone et alJ 120001 : 
lMukhopadhvavll2006l) . it is sensible to consider its effect in 
model studies. This collisional disruption limit, combined 
with the rapid inward drift of m-size particles by which they 
are 'lost' from the local region, led to the concept of a 'm-size 
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barrier' or bottleneck to growth; once large particles exceed 
this barrier, relative velocities become lower, allowing them 
to grow further and to drift less rapidly out of the accreted 
region. Our studies were initially motivated by a desire to see 
if growth into loose fractal clusters and subsequent packing 
could allow the m-size barrier to be crossed. 

In the conceptually simples t models, growth occurs by 
simple sticking of particles dWeidenschillingl 1 19971 120001 
2004; Bluml l2004T) . While it goes against our earthbound 
intuition that macroscopic particles can stick to each other, 
some microgravity and earthbased experiments show that, 
while bouncing transpires at intermediate velocities, stick- 
ing prevails for both low (less than a m eter/second) and high 
(for 13 - 25 m s" 1 ) relat i ve velocities ( Blum & Wurmll2000b 
Blum & Schriipler 2004; jWurm et al.1 120051; iMarshall et al. 
2005; Dominik et al. 2007). Other experiments indicate that 
certain solids (water and methan ol ice, organic mate rial) are 
'stickier' than others (silicates) (Bridg eTet alJI 1996b . How- 
ever, no significant amount of these especially sticky materi- 
als has been found in primitive meteorites. Still, entire chon- 
drites (and by inference entire parent bodies) are composed of 
small silicate objects that seem t o have been gently assembled 
and compacted, at least initially (jMetzler et alii 19921 : iBrearleyl 
1996; Cu zzi & Weidenschilling| l2006): how did this happen? 
The meteorite record (discussed in more detail below) shows 
that many mm-sized solid objects are encased in rims of mi- 
cron and submicron-sized mineral grains. One obvious possi- 
bility is that these rims form by nebula accretion of grains onto 
the underlying core particle s (Nagahara 1984; Metzl eret al] 
199 1 iPaque & Cuzzilfl997L iHua et alJl200l IZega & Busecl 
2003s). In §j2] we show that this dust accretion may be expected 
to occur fractally, leading to porous structures. Collisions eas- 
ily crush this structure, in the process dissipa ting kinetic en- 
ergy a nd allowing colliding particles to stick dBlum & Wurml 
l2008h . 

Recent reviews of the releva nt prop erties of chondrites are 
pro vided bv IBrearlev & Jonesl (119981) . IScott & Krotl d2005al) 
and IWeisberg et alJ d2006l) . Chondrites are dominated by 
mm-size d silicate chondrules, which were melted in the 
nebula dBrearlev & Jonesl 11998b Uones et alJ 120001 120051: 
Lauretta & McSween 200^), but are found in meteorites to 
be embedded in a fine-grained matrix. Formation of chon- 
drules (and chondrite s ) occurred over a p eriod of several Myr 
dScott & Krotl 120053 iKita et al.l 12001 iRussell et al.l [2006k 
iKleine et alJ 120061) . Chondrites can be divided into three 
broad classes - ordinary, carbonaceous and enstatite - with 
each class being further subdivided into more than a dozen 
groups reflecting chemical, mineralogical and isotopic dif- 
ferences. For example, CM carbonaceous chondrites contain 
abundant matrix that has been affected by aqueous alteration, 
while Ordinary Chondrites contain very little matrix and have 
generally incurred only limited aqueous alteration. Violent 
collisional processes occurred after primary accretion which 
affected the contents and appearance of most meteorites, and 
to best understand the primary accretion process one must 
look back through this stage where possible to the rare, un- 
brecciat ed subset of rocks and rock fragmen ts called 'primary 
texture' dMetzler et al Jl 1 992t [Brearleyll 1 993h . 

The dust in chondrites is found to have two physi- 
cally defined components: rims and int e r-chondru l e ma- 
trix dHuss etaD pHt IScottetalJ [19881 IBrearlev! [19961: 
IBrearlev & Jonesl 119981) . Fine-grained rims are clearly as- 
sociated with individual chondrules and other macroscopic 
particles in microscopic images, and usually even stay at- 



tached to the chondrules when the y are disaggregated from 
the host rock dPaque & Cu zzi 1997). Some studies report that 
the composition of these fine-grained rims is uniform across 
a wide range of underlying mineral types, including more re- 
fracto ry (higher-temperature) oxide s which formed much ear- 
lier dBrearleyll 1993b IHua et aLll2002l) and some find dramatic 
variati ons between the composition of rims o n adjacent chon- 
drules dTavloretalJI 1983b IScottetalJ[l984l) . Generally, the 
chondr ules cooled com pletely before accreting these fine dust 
grains (Brearley| |1993l) . Interchondrule matrix, more gener- 
ally dispersed between all the macroscopic components of the 
rock, is also made of fine-grained material. The grain sizes 
in fine-grained rims are noticeably smaller than the ubiqui- 
tously enveloping matrix, even though the comp ositions of 
the rims and matrix are very similar or identica l dAshworthl 
1 1 977b iBrearleyll 1 993b [1991 IZolensky etaD 1 9931) . It has been 
reported that the rim porosity is also s maller than that of the 
surrounding matrix dAshworthl 1 1977b iTrigo-Rodrigu ezet alJ 
2006). The relative abundance of rim and matrix material 
to cho ndrule material v aries from one meteorite class to an- 
other dScott et alj|1988l) : however, the rim mass (or thickness) 
is often found to be pro portional to the mass (or radius) of 
the un derlying chondrule (Me tzler et alJl992blPaque & Cuzzil 
[19971) . 

Several different model efforts have attempted to explain 
some of these properties in the context of nebula sweep-up or 
accretion of the fine -grained rims by chondrules and their like. 
iMorfill et al.l 0998) hypothesized that if a particle had a speed 
relative to the gas which was proportional to its radius, and if 
the chondrules in a region sweep up all the local dust in a one- 
stage event (no ongoing replenishment of dust ), the observed 
rim-core correlation would be obtained. ICuzzil d2004h showed 
that chondrule-size particles in turbulence plausibly exhibit 
just the appropriate (near-linear) dependence of relative ve- 
locity on size, even thou gh most pa rticles obey a square-root 
dependence on radius. Cuzzi ( 2004 ) rel i ed on collisional out- 
comes proposed by Dominik & Tielens ( 1997) for porous ag- 
gregates of fine grains, and suggested that for particles much 
larger than chondrules, the velocity relative to the gas in- 
creases to a point where they ent er an erosion al regim e. 

On the other h and, ISears et alJ d 19931) and 
ITrigo- Rodrigue zet al.l d2006l) question whether fine- 
grained rims are n ebula accretion products at all. 



ITrigo- Rodrigue zet al.l ([2006) point out in particular that 
the fine-grained r ims in CM chondrites, such as discussed 
bv iMetzler et alJ dl992l) . have a porosity of 10-20%, much 
lower than the high -po rosity structures formed by, e.g., 
iBlum & Wurml d2000l) or lBlum & Schraplerl d2004l) . Less is 
known quantitatively about the p orosity of fine-gr ained rims 
in other chondrite types, although Ash worthl (1 1977b states that 
rim porosities are less than 6-15% in ordinary cho ndrites. The 
alterna te that ITrigo-Rodrigue zet alJ ([2P06) and ISears et ail 
(1993) prefer, while they differ in the details, is that the 
fine-grained rims seen in CM chondrites, in particular, are 
created on the parent body from a generic enveloping matrix, 
by some combination of compaction and pervasive aqueous 
alteration. This suggestion might make it harder to explain 
why the grain size is smaller than in the nearby enveloping 
matrix. Nevertheless, the discussion shows that the porosity 
of fine-grain rims is an important diagnostic of their origin. 

In this paper, we develop a detailed collision model 
to study the rimming and accretion processes of chon- 
drules simultaneously and, in a statistical study, quantify 
the growth that can be obtained under a wide range of 
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(uncertain) nebular conditions. Our model treats multi- 
ple components: solid 'chondrules,' submicron grains and 
their very porous nebula aggregates, porous accretion rims 
formed by direct accretion of monomers and aggregates 
onto chondrules, compact rims formed by collisional com- 
pression of pre-existing porous rims, and compound ob- 
jects formed by sticking of rimmed objects, which them- 
selves might become rimmed in dust. Our collisional out- 
comes use physi cal guidelines set by laboratory and the- 
oretical models (iDominik & Tie lensl 119971: iBlu m & Wurm 
I2000t iBlum & Schraplerll2004 lLangkowski et al.ll2007l) . We 
use quantitatively correct closed form relative velocity ex- 
pressions for particle s in turbulence of varying intensity 
dOrmel & Cuzzil [2007b . which capture the increase in rela- 
tive velocity as particles grow by accretion of other parti- 
cles. We use a Monte Carlo approach to calculate the prob- 
ability of different outcomes, over a wide range of nebula 
parameters (level of turbulence, gas and solid density). We 
assess i) the extent to which fine-grained rims can dissipate 
collisional energy and allow growth by sticking to proceed, 
and ii) the extent to which these dissipative collisions com- 
pact initially porous dust rims into lower porosity states. We 
leave for future study the physics of disruptive collisions and 
the details of vertically varying particle density and turbu- 
lent intensity, such as might occur if the global turbulent in- 
tensity is vanishingly sm all (Cuzzietal. 1993; S ekival [19981 : 
iDobrovolskis et alll999tlWeidenschimng|[T997l) . 

This paper is organized as follows. In $2] the collision 
model is discussed. Here we outline the three distinct compo- 
nents with which we model the compound objects that result 
out of accretion of dust on chondrule surfaces and collisions 
with other chondrules. We also introduce the different sources 
of relative velocities particles can obtain in the nebula, and 
calculate the timescales involved in the various accretion pro- 
cesses. We end this discussion with a brief summary of the 
envisioned collisional scenario. §[3]briefly reviews the Monte 
Carlo code with which the coagulation is solved, ^[presents 
the results of our work. First, a few individual models are 
addressed in detail, before we present the results of a param- 
eter study in which many uncertain (mostly nebula-related) 
parameters are varied. In we discuss the effects of a parti- 
cle dominated environment caused by settling of compounds 
on the growth of compounds. We also discuss several obser- 
vational implications, emphasizing in particular the relation 
between the dust in our model to the fine-grained rims seen 
around chondrules in meteorites. We summarize our results 
ind 

2. MODEL 
2.1. Outline 

The central theme of this paper is to model the process of 
dust accretion onto chondrule surfaces and explore whether 
compaction of this dust during inter-chondrule collisions acts 
as a sticking agent, with which significant growth can be 
achieved. Fo r dust aggregates this compaction mechanism is 
well known (IDominik : & Tielensl[l997t [Blum & Wurmll2000t 
IWada et al.ll2007l) : by restructuring of the constituent grains 
the excess collisional energy is dissipated. For dust-rimmed 
chondrules we argue the situation is analogous, except that 
part of the aggregate's interior is now replaced by a chon- 
drule. The prerequisite for such a scenario is the presence 
of a reservoir of dust that is accreted fractally by the chon- 
drules, prese r ving i ts fluffy structure. Here, we follow the 
Mor fill et all d 19981) 'closed box' scenario in which a fixed 



TABLE 1 
List of frequently used symbols 



Symbol 



Description 



n 

A 

Vm.VT 

<P,<P?CA,<i 

<X 

Tf 

PX 

Pd,Pg 

P?.pf 
PUP2 
ij 



Av or A»i2 relative velocity 
fi local Keplerian orbital freque ncy 

a turbulent strength parameter ( 32.2.1) 

y surface energy density 

S fractal growth parameter ( 32.4. It 

e size ratio (e < 1 ) 

nebula pressure param eter (e q. (JJ) 
mean free path (gas) ( 32.2. 1) 
molecular/turbulent kinematic viscosity ( 32.2. 1) 
filling-factor (PCA/porous dust) (Tabled 
collisional cross section 
friction time (eqs. (3] andjSJ) 
gas density over MSN ( 32.2. 1) 
spatial dust/gas density ( 32.2.1) 
specific material chondrule/dust density (Tablef5J 
internal particle density 
dj collision rate betwe en part icles i, j (5]3) 

E collisional energy ( 32.3.21 

£Voll,£br monomer rolling /breaku p energy (eq. 1121 ) 
H g gas scaleheight ( 32.2.4) 

R heliocentric radius 

Re Reynolds number ( 32.2. 1) 

Kg C spatial gas-chondrule mass-ratio ( 32.2. U 

spatial chondrule-dust m ass rati o ( 32.2. It 
St particle Stokes number ( 32.2.3) 

P pressure 
T temperature 

V geometrical (total) compound volume (e q. 1201 ) 

Kg, Opg gas and particle-gas turbule nt velo city ( 32.2.4) 

a geometrical (total) radius ( 32.2.3) 

a p reduced radius, aiaiKf ii + ai ) 

a<i ust dust aggregate radius ( 32.4.1) 

f/ch chondrule radius 

ciq monomer radius ( 3T2~T1 

c g sound speed (gas) ( 32.2.1) 

/comp required dust fraction at collision for sticking (eq. 1161 ) 

fi compound dust fraction by mass (TablefJ) 

/ gco geometry factor ( 32.3.3) 

f p compound porous dust fraction with respect to total dust mass 

(TablefD 

m, md,H! c h (dust/chondrule) mass 

mw-.. or (..),„ mass-weighted averages (see eq. 1301 ) 

? dd> ? dc. f cc dust-du st/du st-chondrule and chondrule-chondrule collision 

times ( 3T4l 

fL / t s large/small eddy turn-over time ( 32.2.4) 

«K orbital (Kepler) velocity 

» r particle radial drift velocity (eq. (8)) 

Note. — List of frequently used symbols that re-occur over various sections in 
the paper. 

amount of dust is injected instantaneously to the chondrule 
population. The compound objects (or, simply, compounds) 
thus obtained are modeled to consist of three phases: chon- 
drules, compact (i.e., restructured) dust, and porous (i.e., frac- 
tally accreted) dust. The restructuring mechanism also holds 
for collisions between compounds, again at the expense of 
the porous phase. In this way a coagulation process is initi- 
ated by which chondrules are accreted into large compounds. 
This coagulation is only stopped when the compounds run out 
of porous dust, such that the collisional energy can no longer 
be dissipated; at this stage, all the free-floating dust has been 
accreted and no more porous dust remains. 

The environment in which these processes take place is a 
key ingredient that enters the coagulation model. A violent, 
low dust density environment leads to high velocity collisions 
which quickly pack down the porous rim, limiting its capabil- 
ity to allow further sticking, while even higher velocities will 
lead to break-up of compounds. On the other hand, if (rela- 
tive) velocities are modest and remain so during the phase in 
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TABLE 2 

Gas and dust parameters 



(1) 


(2) 


(3) 


(4) 


gas density* 


Pg 


28/2.4/0.16 


10 11 g cm 3 


sound speed 


c g 


10./7.6/5.6 


10 4 cms" 1 


mean free path (gas) 


A 


6.9/82/1230 


cm 


temperature 


T 


280/162/89 


K 


pressure parameter 


V 


1.6/2.7/4.9 


xl0~ 3 


large eddy turn-over time 


t L = fr 1 


0.16/0.82/5.0 


yr 


turbulent strength parameter* 


a 


io- 4 




inner eddy turn-over time* 


U 


1.3/12/131 


10 3 s 


gas-chondrale ratio* 




100 




gas-dust ratio* 


Red 


100 




chondrule mean size 




300 




monomer dust size 


"0 


1 




fractal growth parameter 


6 


0.95 




surface energy density dust* 


7 


19 


erg cirT 2 



Note. — Parameters characterizing the gas (upper rows) and 
dust/chondrules (bottom rows). Gas parameters correspond to a minimum 
mass solar nebula (MSN) model of total gas mass of 2.5 X 10~ 2 M Q inside 
10 AU with power law exponents of, respectively, —1,0 and -0.5, for the sur- 
face d ensity and temperature structure as function of radius I Takeuchi & Lin 
I2002D . Columns denote: (1) parameter description; (2) symbol, 7? gc = p g /p c 
and = pdpd', (3) corresponding value with multiple values denoting con- 
ditions at 1, 3 and 10 AU, respectively (3 AU is the default); (4) unit. Param- 
eters indicated by a ★ are variables (default model values are given). 

which compounds accrete other compounds and porous dust, 
this could lead to significant growth. In this work we use vari- 
ous sources for relative velocities: thermal, turbulent and sys- 
tematic. The relative velocity is further determined by the 
internal structure (density) of the compounds, which affects 
their coupling to the gas. The internal structure of compounds 
is reflected in the definition of the 'geometric size': the size 
that corresponds to the effective aerodynamic cross section. 
The evolution of the internal structure during the accretion 
process is therefore of key importance: i.e., do collisions fol- 
low 'hit-and-stick' behaviour in which growth proceeds frac- 
tally, or do collisions keep the filling factor constant. 

2.2. The turbulent nebula 
2.2. 1. Disk physical structure and model components 

The start of the model is defined as the point at which pop- 
ulations of chondrules and dust interact. Before this point, 
the two populations were either spatially isolated or one did 
not exist. This paper is not concerned with the history of the 
two populations - specifically, we avoid the nagging chon- 
drule formation question - but merely define the zero time of 
the model (f = 0) as the point where the two populations mix. 
Of course, the history of the two populations determines to 
a large extent the conditions that prevail at the start of the 
model, e.g., the size of the dust particles. We take 1 /mi 
as the monomer size (i.e., the smallest constituent size of a 
dust grain), which roughly cor responds to the s izes of the fine 
grains observed in chondrites (Ashworthl l 1977b . Thes e grains 
may h ave formed by condensation onto seed grains dChiangl 
2004) and could have aggregated into larger (fluffy) dust par- 
ticles before they interact with the chondrules. For instance, if 
dust condenses out in a region of (recently formed) chon drules 
(IWasson & Trigo-Rodrigue3l2004t IScott & KroJi2005bl) . the 
size distribution will be dominated by monomers. However, 
if chondrules encounter a dust cloud only after a certain time 
since its formation (~ 10 2 ' 3 yr), larger aggregates are ex- 
pected to have formed through monomer collisions. In this 
work, though, the size of the dust aggregates (adust) is simply 
treated as a free model parameter ( §12.41 i. 



Both solid components find themselves in the gas- 
dominated protoplanetary d isk. We use a minimum mass 
solar nebula (MSN) model (Hava shHll981h to determi ne the 
gas parameters (Table [2] after iTakeuchi & Linl |2002|) . At 
R — 3 AU this corresponds to a spatial density of p g = 2.4 x 
10~ n g cirT 3 , a thermal speed of c g = 7.6 x 10 4 cm s _1 and 
a mean free path of A = 82 cm. Studies of chondrule forma- 
tion i ndicate, however, that gas densities may be much hi gher 
(e.g jDesch & Connolly|2002UCuzzi & Alexanderll2006l) and 
we will therefore treat p g as a free parameter, denoting with 
px the density enhancement over MSN. We further assume 
the gas surface density (Z) and temperature (T) profiles are 
power-laws of heliocentric radius (R) and fix the exponents at 
— 1.0 and -0.5, respectively. This choice is consistent with a 
steady disk in which 1 = and an accretion flow, M, inde- 
pendent of radius. 1 The gas-to-chondrule and the chondrule- 
to-dust density ratios, 

^ g c = p g /p c ; and 7? c d = p c /pd, (1) 

are free parameters; for example, 7? c d = 1 means that chon- 
drules and dust are pres ent in equal prop ortion (as observed in 
some meteorite groups, IHuss et al.l2005T) and ft gQ = 100 is the 
standard gas to solids ratio. In our model chondrules follow a 
log-normal size distribution, for which we take parameters of 
300 jt/m (th e mean) and a (log-normal) width paramete r of 0.5 
(see, e.g.. Hubin & Keillfl984l: iNelson & Rubinl 12001) . This 
makes the mean, mass-weighted size of the chondrule popu- 
lation, (a) m = 720 pirn. Note that the initial distribution is not 
necessarily equivalent to the distribution that ends up in me- 
teorites, or the distri bution extracted from meteorites by the 
thin section method (Eisenhou rll996l) . An overview of all pa- 
rameters characterizing the gas and solids is given in Table|2] 
We assume that the gas in the disk is in a turbulent state of 
motion. After Shakura & Sunvaevl ( 1 1973b . the turbulent vis- 
cosity is parameterized as 

vt = ac g H„ = ac 2 JQ., (2) 

where H g is the scaleheight of the gas disk, Q the local (Keple- 
rian) rotation velocity, and a a sca le parameter that determine s 
the strength of the turbulence (Shakura & Sun vaevl Il973h . 
Values for a are very uncertain. If the magneto-rotational 
instability is active it may be up to 10~ 3 dBalbus & Hawleyl 
1199 It iHawley & Balbus|[l99lh : however in regions of low 
ioniza tion it can be much lower (Gammi lEHl [sinoetaD 
2000). The extent of the turbulence is determined by the 
Reynolds number, Re, defined a s Re = Fr/Vm, w ith v m the 
molecular viscosity, v m = c g A/2 (ICuzzi et al.ll 1 993b . The tur- 
bulent spectrum consists of eddies characterized by a scale 
({), velocity (u), and turn-over time (r), between an outer (or 
integral) scale L and an inner (or Kolmogorov) scale £ s . Fol- 
lowing previous works, t^, the largest eddy turn-over time, 
is taken equal to the inverse orbital frequency, ?l = 
and v, = a 1/2 c> (e.g. [Pubrulle et alJ [l995l: ICuzzi et a"Dl2001b 
ISchrapler & Henningll2004l) . The eddy properties at the tur- 
bulence inner scale then follow from the Reynolds number: 

t s = Re~ I/2 f L ; £ s = Re~ 4/3 L; v s = Re^V (3) 

2.2.2. Thermal motions 

When gas molecules collide with a larger (dust) particle 
momentum is transferred, changing the motion of the dust 

1 Provided the a-turbulence model (e g. 0) is assumed and boundary con- 
ditions are neglected (see, e.g., Pringle 1981). 
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particle. These kicks occur stochastically, resulting in a ve- 
locity behaviour known as Brownian motion. The ensuing 
velocity difference between two particles of mass m\ and 
is highest for low masses and high temperatures, 



8k%T(m\ + mi) 



(4) 



where is Boltzmann's constant. For micron-sized parti- 
cles Brownian velocities are a few mm/s; but since Adrm de- 
creases with the -3/2 power of the size of the smallest parti- 
cle, it quickly becomes negligible for larger particles. 

2.2.3. Systematic motions 

The key parameter that determines the coupling of solids to 
the gas is the friction time, Tf. In the Epstein regime the size 
of the particle, a, is small with respect to the mean-free-path 
of gas molecules, A, and the friction time is given by 



Tf 



4c g Pg 7Tfl 2 



(a < U) 



(5) 



For solid 1 /^m-grains the friction time is Tf ~ l0 2 /p X S 
for the default nebula parameters of 3 AU (Table 0, while for 
an a ~ 300 /mi chondrule it takes ~ 10 /px hours before the 
traces of its initial motion are 'erased.' Note that equation (O 
defines a as the geometrical radius of the particle, i.e., the 
radius corresponding to the angularly-averaged projected sur- 
face area of the particle. If the particle is a fluffy aggregate its 
friction time is therefore much less than a compact-equivalent 
with the same mass. If significant growth takes place, parti- 
cles will no longer obey the Epstein drag law; friction times 



are then enhanced with respect to r^ p , i.e 



4a 1 



St _ t" j Bp 

f 9^Re p C f ' 



(a > U) 



(6) 



Here, Re p = 2av pg /v m is the particle Reynolds number, which 
determines the constant C, 2 and v pg the particle-gas velocity. 
Within the physical conditions of the simulations in $4] the 
particle Reynolds number stays below Re p = 1, for which 
C = Re" 1 (I Weidenschillingl fl 9771) . Friction times are then 
independent of v pg . 

One of the well-known problems in the planet-formation 
field is the strong inward radial drift particles of a specific 
size experience, e.g., meter-sized particles at ~AU ra dii or 
cm-sized at ~ 100 AU radii (see Braueretal. (2007) for a 
recent review). This inward radial drift is caused by the exis- 
tence of gas pressure gradients, resulting in a gas velocity that 
is som ewhat less than Keplerian by a difference of magnitude 
tjv k dWeidenschillindl 1 9771: iNakagawa et al Jl 1 9861) with rj the 
dimensionless pressure parameter, defined as 



T) = 



1 1 dP 
'2Rn 2 p~ g dR 



(7) 



However, particles do not experience this pressure term and 
instead attempt to move at Keplerian velocities, faster than 
the gas. The ensuing drag force removes angular mo- 
mentum from the particle resulting in a radial velocity of 



2 Compared to, e.g., Weidenschilling 1 1977) the definition of C has been 
scaled down by a factor of 24. That factor is already present in equation ^6). 
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Fig. 1. — Comparison of systematic and turbulent velocities as function 
of particle Stokes number. Plotted are radial v eloci ties (eq. (8), solid grey 
curve) and turbulent velocities with e = (eq. 1101 . dashed curves). Note 
the change in slope between the linear and the square-root turbulent velocity 
regime which happens at St = Re -1 ' 2 . Turbulent velocities dominate over 
systematic velocities for St < a/rj (provided this is ^ Re" 1 ' 2 ). The Stokes 
number corresponding to a 300 /.im-sized chondrule at the default nebular 
conditions of 3 AU (see Table |5J is also indicated. All velocities peak at 
St = 1 (light-dashed vertical curve). 



(IW^denschilli i^[T977l) 3 



u, = -- 



2St 
1 +St 



T]Vk, 



(8) 



where we have defined St = Tf Q. This systematic radial drift 
velocity peaks at St = 1 . Chondrule-sized particles, however, 
are generally sufficiently well coupled to the gas (St <K 1) 
so that radial motions are low in most physical conditions; 
turbulent-induced motions then dominate (unless a is really 
low). However, when particles grow in size, systematic mo- 
tions may take over from turbulent velocities (see Fig. |TJ. 

2.2.4. Turbulent motions 

For a Kolmogorov spectrum, turbulence leads to mean 
(large scale) velocity fluctuations of v g = (3/2)'^^ = 
(3/2) 1/2 a 1/2 c g dCuzzi & Hoganll2003l) . Due to their inertia, 
solids do not instantaneously follow these fluctuations but re- 
quire a time Tf before their motions align. This leads to a net 
relative motion, p pg , b etween the gas and the solid particle of 
dCuzzi & H ogan 200l 



St 2 (Re 1/2 - 1) 
(St + l)(StRe 1/2 + 1) 



(9) 



where the Stokes number, St, is the ratio between the friction 
time and the large eddy turn-overtime, i.e., St = Tf/iL = TfQ. 4 
The limiting expressions of equation d5). fpg = Stu g for St <s 



1 and v, 



for St » 1, respectively, correspond to the 



3 In equation {8) we have not accounted for collective effects when the 
particle densi ty is comparable to or higher than the gas density. Equation (8) 
then changes I Nakagawa et al. 1986). Angular momentum exchange between 
the dust and gas dominated layers I Youdin & Chiang 2004J) is an other process 
to be accounted for, but its significance is relatively modest I Br auer et all 
12001) . 

4 St = Tf/fL is the formal definition for the Stokes number. In the a- 
turbulence model ?l = ^ 1 and the definitions for St in equations (8) and j9) 
coincide. 
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cases of particles that are well coupled (small particles) and 
poorly coupled (larger particles) to the gas. 

The calculation of particle-particle relative velocities does 
not follow directly from the y pg 's since particle velocities can 
become very incoherent in turbulence (e.g., Ai>i2 + |i>i g - U2gD- 
Consider, for example, two small particles entrained in the 
same eddy. If their motions are aligned, no relative ve- 
locity is present; it is only within a time Tf after being 
caught in the eddy that these particles have the chance to 
develop relative motions, provided their friction times dif- 
fer (i.e., T\ + T2). The problem of finding suitable (i.e., 
closed-form) expressions for Au\2 is importa nt since these are 
key to any model of du st coa gulation (e.g.. |Weidenschilling| 
19971; ISuttner & Yorkell20qit iDullemond & Dominikl 120051: 
Nomura & Nakagawal 120061) inclu d ing th is work. Follow- 
ing earlier works of IVolket al.1 (1 19801). I Weidenschilling 
dl984.lMarkiewicz et all d 199 lb and lCuzzi & Hoganl d2003l) . 
Ormel & Cuzzil (120071) have presented closed-form analytical 
expressions for Aun (with a margin of error of ~ 10%) in 
terms of the Stokes number of the particles: 

2 fRe 1/2 (Sti-St 2 ) 2 forn^ 
= \[2yt-(l + s) + t L ( 1 4 + ;g s )]St 1 forS^r, 

(10) 

In these expressions t\ (or Sti) always corresponds to the par- 
ticle of the largest friction time and s = tiIt\ < 1. Near 
the T; = f s turning p oint the expression for Avn is somewhat 
more complex (see Ormel & Cuzzi 2007]). j/* is a numeri- 
cal constant of value y* d ^ 1.60 if t\ <k t\_.\ however, when 
Ti ^ /l h becomes a function of t\ and drops to unity at 
ti = ?l- In that case we approximate y\ by an interpolation 
function. In Fig.Q]Aui2 is plotted for two values of a in the 
limit of St = St; » St2 (dashed curves). The three regimes 
of Equation ( fTOb are clearly distinguishable: the linear regime 
for Tf < f s (or St < Re~'^ 2 ; small particles); the square-root 
regime, f s < Tf < t\_; and the high Stokes regime, St > 1, 
where particles decouple from the gas. Small particles like 
chondrules have Stokes number St « 1; whether they fall 
into the linear or square-root velocity regime depends on their 
sizes in relation to the gas parameters (e.g., a,p g ) of the disk. 

Figure Q] also shows the systematic drift velocity (eq. JSJ, 
solid curve). If one assumes that St2 <sc Sti the radial drift 
curve also gives the relative velocity a particle with St = Sti 
has with a much smaller particle. Actually, for Sti < 1 turbu- 
lent and systematic drift relative velocities (eqs. [8) and IfTOl ) 
depend only weakly on St2 (the lower Stokes number); the 
curves in Fig.Q]can therefore be interpreted as the typical rel- 
ative velocity a particle of a given Stokes number has with 
particles of similar or lower Stokes numbers. Figure Q] shows 
that for very small particles (Tf <K f s or St <K Re~'^ 2 ) turbu- 
lent velocities only dominate when a > 10 . Then, when 
turbulent velocities flatten out in the square-root regime, the 
radial drift motion may catch up with the point of intersection 
lying at St ^ a/rj, provided a is not either too low or too high. 
For any model with a < 77 radial drift motions will eventually 
dominate: a regime of high relative velocities (~ 10 m s~') is 
therefore unavoidable. 

2.3. Collisions between dust-rimmed chondrules 

Theoretical studies and laboratory experiments have shown 
that the outcome of grain-grain or aggregate-aggregate col- 
lisions dep ends on its ratio of the kinetic energy to a criti - 
cal energy (Domini k & Tielens|[l997t iBlum & Wurmll2000l) . 



TABLE 3 

Material properties of compounds phases 



Phase 


Specific density 


Mass fraction 11 


Filling factor 


chondrule 


p ( c s) = 3 g cm- 3 


(l-/d) 


<f>ch = 1 


compact dust 


(S) n -3 

P\ = 3 g cm 


/d(l-/p) 


= 0.33 


porous dust 


(s) 1 -3 
P' ( = 3 g cm 


fdfp 


P d < 0.15 



a with respect to entire compound (sum equals 1); /j = mass fraction 
in dust; / p = porous mass fraction of the dust. 




chondrule compact dust porous dust 



Fig. 2. — (A) Schematic representation of the three-phase model of dust- 
rimmed chondrules: chondrule (dark), compact dust (grey) and porous dust 
(pattern). The inset (B) shows the substructure of the dust that consists of 
micron-sized monomers. (C) If chondrules collide, the collision energy is 
dissipated by transferring dust from the porous to the compact dust phase. 
Figure not to scale. 

Specifically, porous dust is accreted when chondrules sweep 
up dust grains, or aggregates of dust grains, at collisional en- 
ergies (E) that stay below the energy threshold for restruc- 
turing, 5E ro \i where E m u is the energy required to roll one 
contact area over the surface of the grain. This can lead to 
a very open structure of filling factors (0 p d) that are lower 
than the filling factors obtained in particle-cluster aggrega- 
tion, 4>pca — 0.15 (see below, {32.3.11 ). As outlined in 32.11 
the absence of restructuring during dust sweep-up is crucial 
since the resulting porous structure can then be compacted in 
more energetic collisions - the collision of compounds - pro- 
moting further growth. In compound collisions, E > E m u 
and the dust within the compound will restructure, dissipat- 
ing a unit of ~ E m u for each dust grain that is involved in the 
rolling motion. The porous dust that is involved in restructur- 
ing compacts to a higher filling factor, </> c( j ( 32.3.11 ). 

Compounds can then be represented as a three phase struc- 
ture: chondrules, compact dust and porous dust. Two num- 
bers, the dust fraction f& and the porous dust fraction f p , quan- 
tify the relative importance of each phase within a compound 
(see TablefJ]). The internal structure of each phase is further 
characterized by its filling factor, 0. A schematic picture of 
the structure is given in Fig. [2] 

2.3.1. Collisional compaction 

The accreted dust mantles surrounding chondrules can have 
a very porous and fractal structure. Typically, grains in these 
rims will be bonded to two other grains in large string-like 
structures. When two rimmed chondrules collide, contact will 
be established between two (or a limited number of pairs of) 
grains and these grains will bear the full brunt of the collision. 
Once the force on these grains exceeds the critical rolling 
force, they start to roll (restructuring). The rolling of these 
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Fig. 3. — (solid curve) Compression of PCA aggregates, obtained by 
static compression of 'dust cakes' created by random ballistic depositi on of 
a = 0.75 /an Si0 2 spheres JBlum & S chrapler 2004; Blum et al. 20®. The 
uncertainty in the measurements is denoted by the grey area. Values for for 
several amounts of packing configurations are shown. (Data on the compres- 
sion curve experiments kindly provided by Jiirgen Blum.) 

grains may enable contact formation between more pairs of 
grains, thereby promoting compaction and at the same time 
reducing the force per contact. Compaction will stop when 
the force on newly made contacts drops below the rolling 
force. Compaction may also stop because the resulting struc- 
ture is too rigid to allow for further rolling, i.e., the rolling 
grain made contact with too many grains. Since forces are 
propagated through such compacted structures, this means 
that none of the grains involved experiences a force exceed- 
ing the rolling force. The compression of the contact area in 
a collision between two monomers will give rise to an elas- 
tic repulsion force slowing down and eventually reversing the 
collision. The absolute value of the repulsive force will be 
set by the kinetic energy of the collision; in the Hertzian limit 
the sum of the forces on the individual contacts scales with 
the remaining kinetic energy to the 3/5th power. For collid- 
ing aggregates, rolling of the contacting monomers provides 
an additional energy dissipation channel. But once the result- 
ing structure is too rigid to enable further rolling, compression 
of the individual contact areas will provide the repulsion re- 
quired to slow down and possibly even reverse the collision if 
the collision is energetic enough. The collision partners then 

b ounce. 

Blum & Schrapler (2004) have designed experiments to 
measure the compaction of dust cakes under uniaxial com- 
pression. These dust cakes were grown through a Particle- 
Cluster- Aggregation (PCA) method by deposition of individ- 
ual 1 .5 yum diameter monomers at low velocities where rolling 
is not a factor and growth occurs through a hit-and-stick pro- 
cess. The volume filling factor of the resulting aggregates was 
measured to be 0. 15, in good agreement with numerical simu - 
lations of this process (Ossenkopl [l99lrvVatson et al.lll997l) . 
At this volume filling factor, the typical coordination number, 
i.e., the number of neighbours wi th which the monomer i s 
in contact, is calculated to be 2 ((Ridgwav & Tarbucklll967l) . 
These dust cakes were exposed to a unidirectional pressure 
in a static experiment. Figure [3] shows the resulting volume 
filling factor as a function of t he applied unidirectional pres- 
sure dBlum & Schraplerl l2004). The results show that com- 



paction is initiated at an applied pressure of ~ 10 dyn cm . 
If we assume that the number of monomers per unit area being 
pressed on is given by, 

N/A = (nal^y 1 ~ 2 x 10 8 cm" 2 (11) 

with flo the radius of the monomer, the force on an individ- 
ual monomer becomes ~ 5 x 1 0~ 5 dyn. This is very clos e to 
the roUing force of 7 x 10~ 5 dvn lBlum & Schrapler! d2004l) ex- 
trapolated from the measurements by lHeim et al.l (Il999t) . As 
the dust cake compacts and the average coordination number 
increases, the structure becomes more rigid and resistant to 
further compression (under these conditions, see below). Fi- 
nally, at a pressure of ~ 10 6 dyn cirT 2 the structure is dense 
enough for rolling motions to be inhibited. This corresponds 
to an average coordination number of 3.9 and a filling factor 
of 0.33 (Fig. 0. 

The conditions in the protoplanetary disk under which dust 
rims are formed by grain accretion and under which they 
evolve through collisions with other rimmed grains differ 
from those in these laboratory studies. First, the initial struc- 
ture of the dust rims accreted on the chondrule surface may 
differ from PCA. Although low velocity collisions are ex- 
pected, the monomers may collide preferentially among each 
other before colliding with a chondrule. In that case the re- 
sulting structure is referred to as Cluster-Cluster Aggregation 
(CCA), a process that leads to much lower filling factors than 
PCA. Whether PCA or CCA is preferred depends on the rel- 
ative abundance of dust and chondrules and their relative ve- 
locities. So, the filling factor of the porous dust component 
may start lower than the experimental one in Fig. [3] How- 
ever, since this process is directly tied to the rolling force ex- 
perienced by the monomers that make contact, we expect that 
this difference in initial structure will have no influence on the 
critical pressure required for the onset of compaction. We ex- 
pect, likewise, that uniaxial compression of dust rims grown 
by CCA will stall at 0.33 filling factor since this is again a 
property of the resulting structure; e.g., at these kinds of vol- 
ume filling factors, monomers in the dust rims will have been 
organized in 'stabilizing' structures. 

However, under nebular conditions continuous impacts will 
arise from random directions; it is therefore likely that colli- 
sional compaction under these conditions will result in higher 
volume filling factors than the unidirectional compression ex- 
periments would indicate, possibly even as high as 0.5 (the 
value of (p = 0.64 characterizes for Random Clo se Packing, 
RCP) We note that lBlum & Schrapler! d2004l) and lBlum etal] 
(2006) in their compression studies did approach RCP when 
applying an omnidirectional pressure of ~ 10 9 dyn cm -2 . 
Note, however, that omnidirectional pressure is not achieved 
in collisions between two bodies in an open environment; i.e., 
the dust has the chance to spread perpendicular to the di- 
rection of compression, and the obtained high volume fill- 
ing factor may not be generally attainable. Indeed, for 
(f> > 0.33 rolling motion become impeded and we do expect 
that in order to 'crush' the dust (rims) to RCP values much 
higher pressures are required. Studies in dicate pressures of 
~ 10 10 dyn cm" 2 in orde r to reach RCP dMartin et alj|2003b 
iTanwongwan et al"]l2005h . This second stage of compaction 
would correspond to a very different collision regime charac- 
terized by much higher energies. Adopting the Hertzian limit, 
we expect that this higher pressure, c.q., force in the contact 
area corresponds to a impact energy which is higher by a fac- 
tor (10 9 /10 6 ) 5 / 3 or a collision velocity higher by a factor 10 5 ^ 2 
over the velocity/energy required to initiate compaction. In 
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TABLE 4 
Critical energies 



Expression 




Breaking 




Rolling 


Theoretical a 




= A brr 5 V 3 /£ 2/3 


Eroll 


= A TO nfcrit^y 


DT97 Prefactors" 




a 9.6 


Aroll 


= 59 


Measured energies^ 


Ebt 


= 1.3 X 1(T 8 erg 


firoll 


= 1.7 X ltr 8 erg 


Empirical prefactors c 




= 2.8 x10 s 




= 1.8 X 10 3 



Note. — Comparison between predicted and measured critical energies 

for breakup and rolling. 

" Theoretically derived express ions for E bl , E m [[ f rom Chokshi et al] 11 9931) 
(for the breakup energy) and Dominik & Tielens) II 9971) (for rolling) and 
corresponding pre-factors, A^ r , A ro u . We define f cr ; t = 10~ 8 cm. 
b Values adopted from Blum & Wurm ( 2000) for parameters of uq = la^ = 
9.5 X Kr 5 cm, y = 19 erg cm -2 and £ = 3.7 X 10" dyn cm' 2 . The origi- 
nal measure ments were perfor med by Poppe et al. ( 1999) (for the breakup 
energy) and Heim et al. 1 1999) (for the rolling energy). 
c Empirically derived prefactors from the theoretical expressions with the 
measured values for Eb T and E m \\ . 

this study, while acknowledging that higher filling factors are 
plausible, we have for simplicity assumed that <p = 0.33 is the 
limiting value. 

2.3.2. Acquisition of a porous dust layer 

Two critical energies - the breakup and rolling energy - reg- 
ulate the behaviou r of the dust (porous a c cretion/compaction) 
upon collision: (IChokshi et al.l 119931: iDominik & Tielensl 
fl997tlBlum & Wurmll2000l) ~ 



5/3 4/3 



£2/3 



E mU = 67r 2 £ cr ita^y = A^^a^y, 



(12) 
(13) 



where = 0102/(01 + 02) is the reduced radius of the col- 
lision partners, y the surface energy density of the material 
and & Young's elastic modulus (assuming the same materials 
collide). £ cl i t in the £ ro u expr ession is some critic al distance 
used to initiate rolling which Dominik & Helens] (1 1997b as- 
sumed to be on the order of the atomic size, £ cr j t = 10~ 8 cm. 
Using thes e definitions th e con stants Ab r and A m \\ are dimen- 
sionless. iBlum & Wurml (120001) have experimentally deter- 
mined the breakup and rolling energies (see Table B) and 
found these to be higher than the IDominik & Tielensl (1997) 
theoreti cal predictions. How ever, apart from a scale fac- 
tor, the Blum & Wurm (2000) experiments agreed well with 
the IDominik & Tielensl (Q997) model; that is, collisions can 
be separated into the regimes of perfect sticking, restructur- 
ing and fragm entation. We therefore ap ply the mechanism 
put forward by Dom inik & Tielensl ( 1997) but use pre-factors 
(Ab r , A,oii) from the experimental results (last row of Tabled. 
Note that for micron-sized particles the rolling and breakup 
energy are similar. 

When two particles meet, direct sticking occurs if the col- 
lision energy, E, is dissipated at the first point of contact; 
i.e., E < Estick, where £ s tj C k is related to the b reakup en- 
ergy as Estict = 0.22£br (IDominik & Tielenslll997l) . Writing 
E = jm^iAv) 2 with = m\m2l(m\ + mi) the reduced mass 
and Au the relative velocity, the criterion E < E % <ac\ translates 
into a threshold velocity of 



2E« 



V 



= V0-45A br 



0.33A b f (p^f 1 V V /2 < V /2 £~ 1/3 







-5/6 


f ^ ) 


\ao) 






3gcm \ 



-1/2 



y 



5/6 



19ergcm~ 2 / \3.7 x 10" dyn crrr 2 



-1/3 



,(14) 



where we have assumed that like materials meet (i.e., same 
y, fi) , and where the reduced mass has been parameterized as 
nif, = N^mo with wo = 4-7rp d j Oq/3 the mass of the (smallest) 
grain. Thus, = 1/2 for equal-size particles and = 1 for 
very different size particles. Equation ( fl4l i shows that micron- 
sized silicate particles have no problem to stick to each other 
at velocities of ~ 10 cm s~' . This also holds for collisions be- 
tween /mi-sized grains and chondrules since it is the reduced 
size 0^ that enters the equation. However, at higher velocities 
the grains will bounce off. 

In collisions between chondrules (oo ~ 300 |/m) the stick- 
ing velocity falls below ~ cm s~' , lower than the velocities be- 
tween chondrules for most values of a (see Fig . |TJ . Also, for 
chondrules, the assumption of a smooth, spherical surface on 
which the physics behind equation ( fT3l is based breaks down. 
Although surface ro ughness increase t he sticking capabilities 
for/im-sized grains (Pop pe et al.l2000h . the asperities in chon- 
drules are probably too large to favour sticking. However, we 
now expect the previously accreted porous dust layer to act as 
the sticking mechanism through a dynamic restructuring and 
compaction of the constituent grains ( ^2.3.3l l. 

Another important collision is between a chondrule (or 
compound) and a dust grain or dust aggregate. Here, the 
criterion for sticking with out restructuring is E < 5E ro \\ 
(Dominik & Tielens 1997). This translates into a velocity of 
(using the same substitutions as above) 



f.s 



st.aggr 



10E ro 



V 



1.6 x 10 2 cms" 1 AT I/2 



ao 



(s) 



3gcm 3 J 



-1/2 



1/2 



r 



go 
Ifmij 



19 erg g" 



v-1/2 
(15) 



in which now > 1 roughly corresponds to the number of 
grains in the aggregate. Thus, at moderately low velocities 
small aggregates will hit-and-stick, preserving their porous 
structure. This 'hit-and-stick' behaviour will also be referred 
to as 'fractal accretion.' However, when aggregates become 
large or a is high (e.g., a > 10~ 3 and px = 1; Fig. [TJ some 
compaction is likely to occur. However, in this study we have 
ignored this effect (for reasons of computational efficiency) 
and simply assumed that all dust accretion occurs fractally. 
Although invalid in a violent collisional environment, the con- 
sequences of this assumption are marginal as the porous dust 
on the chondrule surface is quickly compacted by colliding 
chondrules in any case ( ^2.3.3l l. 

2.3.3. Collisions between dust-chondrules compounds 

lHatzes et al.l dl991l) studied collisions between cm-sized 
particles and found that sticking forces increased significantly 
when a frosty layer was present. While they attribute this en- 
hanced sticking to interlocking of jagged surface structures, 
this effect probably reflects energy dissipation due to restruc- 
turing. In the case of chondrules, rimmed by a layer of fluffy 
dust, the situation is analogous: the fluffy structure allows the 
collisional energy to be dissipated away. Assuming that each 
monomer (of mass mo, size ao and internal density p^j ) in 
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Fig. 4. — Contours of / cora p (dashed curves) - the fraction of the (com- 
bined) compound mass that must be involved in restruc turi ng to dissipate 
away the collisional energy to stick the compounds (eq. 1161 ) - as function 
of collision velocity (.r-axis) and size ratio (y-axis). Equal internal densities 
are assumed, pi = p2 = 3 g cm" 3 and an = 1 fim The criterion for sticking 
is /comp 2 (/d/p)m ( see text). For low velocities or size-ratio's compaction is 
insignificant. Collisions with / comp > ( f<i fn )m com pact all their dust. When 



! complete destruction occurs (see j|2.3.5t . 



the porous dust layer is capable of absorbing an energy E I0 a, 
E/E m u monomers are needed to dissipate the total collision 
energy. Expressed in terms of mass, a porous mass fraction 
of at least / comp must be available, with / comp the ratio of the 
required mass in porous dust to the total mass of the collision 
partners, 



/c 



m E/E, 



roll 



47T 



comp 



m\ + rri2 



,,,2 



(s) 2 - 

Pi «oT 



3.7 x 10 



-2 



m\m,2 1 

(s) 

Pi 

3 g cm - 



/im 



Av 
10 cm s~ 

y 



(Av) 2 



19 erg cm 



-2 



■(16) 



This equation reveals a few important results. First, / cornp de- 
creases with smaller dust grains (smaller an); although more 
monomers are required to dissipate the same collision energy, 
the total mass of the monomers that restructures is less. Also, 
the dependence on velocity is rather steep; at very low ve- 
locities the amount of compacted material is negligibly low 
(/comp ^1), while visible compaction happens in a quite re- 
stricted velocity band (see Fig. |4j». Another important point is 
the mass-dependence in equation ( fTol l. / comp does not depend 
on the absolute masses of the particles involved, but, through 
the m 2 /mi/«2 factor, rather on the mass-ratio of the collision 
partners. Thus, a collision between particles of very unequal 
size has a lower / comp than equally-sized particles colliding at 
the same velocity and, therefore, a higher probability to stick 
(see Fig. |4j. This is of course due to the reduced mass that 
enters the collision energy. 

2.3.4. Collisional recipe 

We will now quantify how the collisions affect the struc- 
tural parameters of the compound, i.e., the /<j and f p quan- 
tities. Equation ( fl6l > gives the mass fraction of the colli- 
sion partners that must be compacted, which, for sticking, 




Fig. 5. — Projection of two compounds colliding at an impact parameter 
b. The ratio of the shaded region (A) relative to the cross section of each 
compound (ncr ) determines the fraction of the porous dust that can be used 
during the collision, i.e., fgeoj = A/ncq. A is obtained by subtracting the kite 
OiBChC from the two circular sectors spanned up by /?[ and pi (eq- 1171 ). 



must be less than the mass fraction available in porous dust 
(/d/ p ), averaged over the collision partners, i.e., (fif P ) m = 
(mi/d,i/ P ,i + >«2/d,2/p,2)/('«i +m 2 ). If /comp < </d/ P )m enough 
porous dust is present to absorb the collisional energy and the 
two compounds stick. A fraction / comp is then transferred 
from the porous to the compact phase. If / CO mp > (fpfd)m> 
however, restructuring cannot dissipate all the collisional en- 
ergy. As mentioned before, we do not include other en- 
ergy dissipation channels but simply consider all collisions 
in which / comp > (fif p ) m to result in bouncing; / p is then set 
to for both particles. This means fragmentation of com- 
pounds or ero sion of the porous rim are neglected (but see 
below, EE). 

A further restriction to the amount of dust that can be com- 
pacted is made when we account for the geometry of the col- 
lision. Then, only a fraction (/ geo ) of the compound (and of 
its dust) is involved in being compacted and dissipating en- 
ergy. We estimate / geo from the intersection between the par- 
ticles' trajectories. This intersection actually is one between 
a cylinder and a sphere, but here we approximate it as a 2d 
intersection between two circles which meet at an impact pa- 
rameter b. The area of the intersection, A, can be calculated 
by elementary geometry as (see Fig. [5]) 



A(«i, A2. b) — 2a\ arccos 



2ba\ 



+ 2a\ arccos 



(b 2 +a\-a\^ 
2ba2 



- y(—b + a\ + a2)(b + a\ — a2)(b - a\ + a2)(b + a\ + 41J,) 

from which / ge0j ; for each particle is calculated as / ge0j ; = 
A/naf. Equation dTvT > is valid for impact parameters of \a\ - 
flfel < b < a\ + A2- For b <\a\ — a^\ the intersection equals 
the projected area of the smaller particle, while A — for 
b > a\ + CL2- The total mass-fraction of the particles that 
participates is (mi/ d i/pi/ g eo,i + m 2 fi2fp2f ge0 ,2)/( m i + ™ 2 ) and 
this has to be greater than / comp for sticking. Note, how- 
ever, that inclusion of the / geo factor might be too restric- 
tive: since the sound sp eed inside aggregates (~ 30 ms~'; 
Paszun & Dominik 2008) is usually much higher than Ay, the 
energy will be quickly transferred along movable structures. 
For completeness we therefore consider both extremes: / geo 
determined from equation ( TTTb and / geo = 1 . 
When the collision results in sticking, the total dust and 
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porous dust fractions are calculated as follows, 



/d = 



fd\™\ + fdim 2 _ 



fv 



mi + ni2 

fdifpim + fdifpim - /compel + mj) 
fl(m\ + mi) 



(18) 
(19) 



For example, in equation (TT9b the three terms in the numer- 
ator denote, respectively, the mass in porous dust of particle 
1, the porous dust mass of particle 2, and the porous mass 
transferred to the compact dust phase. In order to get the rel- 
ative amount of porous dust this is divided by the new total 
dust mass (calculated in equation ( TT8l ) of the compound. If 
the collision results in a bounce, /j stays the same for both 
particles and / p is reduced by a factor 1 - / geo . 

2.3.5. Role of fragmentation 

As the above formulas for /j and f geo suggest, fragmenta- 
tion is not included in the collision model. The large number 
of particles produced by a fragmenting collision is especially 
problematic in the context of a Monte Carlo program, where 
the num ber of particles is limited ( see @. However, the re- 
sults of Dominik & Tielens (1997) provide insight into the 
stage at which fragmentation becomes important. They found 
that the onset of fragmentation starts at E a 0.3./V c .E'br with 
catastrophic disruption at energies of a lOA^Zibr, where N Q 
is the total number of contacts (roughly equal to the number 
of monomers, N) in an aggregate. Recalling from 32. 3. 21 that 
£br ~ E io h the catastrophic fragmentation limit corresponds to 
/comp ~ 8 in terms of equation ( TToT l. The corresponding curve 
in Fig. |4] then shows that fragmentation becomes important at 
velocities above a few 10 2 cm s _I (for equal sized particles; 
as the mass disparity increases the fragmentation velocity in- 
creases). The m/s transition for the onset of fragmentation 
is in agreement with previous studies (Blum & Miinchl [19931: 
iBlum & Wurml 12000) but compact structures at high filling 
factor may require more effort to fragment than their fluffy 
counterparts (D. Paszun, priv. comm.). From Fig. Q] the crit- 
ical velocity can be translated into a Stokes number. We will 
a posteriori check in which models fragmentation is expected 
to play a dominant role. 

2.4. Evolution of the internal structure 

The aerodynamic properties of the compounds, which de- 
termine their coupling to the gas, alter with accretion of dust 
and mutual collisions between compounds. These properties 
are quantified in the friction time, Tf (see eq. 15)), of the par- 
ticles - essentially a measure of its mass-to-surface area ratio, 
m/na 2 . It is an important parameter since lower friction times 
mean lower relative velocity between the particles ( 92.21 >, and 
low relative velocities between compounds, in turn, imply po- 
tential to grow large. 

Through the adopted three phase model of compounds 
( 92.31 1 the total geometrical volume, Vg eo , defined by a, can 
be reduced to its three components: i) chondrule, ii) compact 
dust, and Hi) porous dust, i.e., 



yfl 3 = V = V ch + Vcd + Vpd- 



(20) 



In 92.3l the filling factors of the chondrule and compact dust 
phases have been fixed at 1 and 0.33, respectively, so that V; 
for these phases linearly corresponds to the mass inside these 
phases. However, for the porous dust phase, this does not 
have to be the case: the dust fluff-balls that are accreted can 



be of different size and porosity. Besides, if the compound 
accretion process itself proceeds fractally, the porous phase 
becomes a mixture of porous dust and voids created by the 
hit-and-stick packing of the compounds. This could lead to a 
much reduced filling factor of the porous phase (see 92.4.31 1. 

With these issues in mind we discuss the three accretion 
modes that are at work and outline the implications for the in- 
ternal structure of dust and compounds. In a largely chrono- 
logical order these are: i) dust-dust, ii) chondrule-dust, and 
Hi) compound-compound accretion. 

2.4. 1 . Dust-dust aggregation 

If the dust distribution initially consists of monomers of size 
flo, the number density of dust particles is likely to be much 
larger than that of chondrules. Therefore, monomers proba- 
bly form aggregates before they are themselves accreted by 
chondrules or compounds. Provided the collisional energies 
involved stay below 5£ ro ii (see 92.3.2b . the dust aggregates 
will hit-and-stick, leading to fractal growth (Meakin & Donnl 
Il988t IBlum et aUfeOOOl) . In addition, if this collision is be- 
tween similar-size aggregates, the outcome is referred to as 
cluster-cluster aggregation (CCA). In that case fractal aggre- 
gates form with a surface area that grows faster with mass 
than in th e compact case, i.e, ggj^ « m 6 , or, V oc m 35 ^ 2 with 



6 = 0.95 dOssenkopfl 1 993t iPaszun & Dominikl2006l) as com- 
pared to V oc m 1 ' 3 (or 6 = 2/3) for compact growth. Using the 
relation V = VoN 36 ^ 2 , where Vo denotes the volume of a sin- 
gle monomer and N the number of monomers the aggregate 
contains, the filling factor evolves as 



NVp 
V 



2/3<5-l 



/ \ 2/5-3 
/ Adust \ 

\ A() / 



N 



1-35/2 



(21) 



We consider two mechanisms through which dust aggre- 
gates can form: i) Brownian motion and ii) differential veloc- 
ities due to turbulence. For simplicity, equal particle sizes are 
assumed at all times. The timescales involved are determined 
by the particle number density (n d ), size (adust) an d relative 
velocities (Au) between the particles, i.e., 



fdd = (n d o-Avy 



1 n W 

1 Pd fl dust« 



(22) 



4 ™L 



3 Pd An 

in which we have used rid = Pd/m, cr = n(a\ + ai) 2 

for the collisional cross section, and m — 47ra 3 |ust 0p[ | SJ '/3 for the 

mass of a dust aggregate, with py the specific material density 
of the dust. The relative velocities in the case of Brownian 
motion and turbulence read, respectively, 



Au BM = 



8k B T 



3k B T 



^ 1/2 ; 



Av l 



"s Pd 



fldust^, 



(23) 



(24) 



Is CgPg 

If Tf < f s , turbulent relative velocities are in fact determined 
by the spread within their friction times (see 92.2.41 ) and a 
numerical factor of, e.g., 0.1 in front of equation (l24l i can be 
inserted if the particle distribution is monodisperse. This ex- 
pression is further reduced by writing vjt s 
For Brownian motion, the timescale fdd then becomes 

n(pf? l2 K %i ^ 3/2fl5/2 ^ L4xl()2yrp -y/2| 



JM 
'dd 



6p g ^3k B T 



Re 1/4 ar 1/2 c g a 



v5/2 
"dust \ 



1 pm ) 



(25) 
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and for turbulence, 
1 



lOx-ftcdRe 

3 s 



4/ V 1/2 £T' * 6xl0 2 yrp. 



\io- 4 ) 



-3/4 



n ch na 2 ch Av 3p c Av 3Re 1/4 a 1 /2Q 



(26) 

where the expressions are evaluated for the default parameters 
of the R = 3 AU model (Table H with 7? gd = ft^Kd = 100 
the gas to dust mass ratio, and where the factor of 10 in 
equation (|26| | follows from the considerations given above. 
For turbulence the increase in geometrical area due to the frac- 
tal growth is exactly cancelled by the decreased Tf so that the 
timescales for turbulence become independent of size and <p 
and growth progresses exponentially; for Brownian motion 
the gr owth of aggregates (in mass) is proportional to t 2 dBluml 
12001 . 

Equations d25l l and d26l l show that aggregate formation is 
initiated by Brownian motion. Turbulence can take over at 
high a but the dust is then also quickly swept up by chon- 
drules. At high gas densities aggregates can grow large. 

2.4.2. Dust-chondrule/compound accretion 

The size of the dust aggregates at which they are accreted by 
chondrules/compounds (adust) depends on the dust-chondrule 
accretion time fdc in relation to fdd- The timescale, fdc, for a 
dust particle to encounter a chondrule of size a c h, is 

1 4p^ ) a ch 47? gc _ m ( a \ 

fdc = — = t— rr = — - — * 240yrp x ^— j 

(27) 

where a monodisperse distribution of chondrules that domi- 
nates the cross section (a c h » adust) and the velocity field is 
assumed, and 7?g C = 100. Equation d27| i again assumes the 
relation Au ~ v^j^/ts, although chondrules might also fall in 
the square -root part of the velocity regime ( {32.2.41 i. Equations 
d26l ) and ( f27b show that in turbulence the dust is preferentially 
swept up by chondrules. On the other hand, in Brownian mo- 
tion, velocities are always determined by the smallest particle 
(the dust); the grains therefore simply collide with particles 
that dominate the total cross-section. Thus, if turbulence dom- 
inates the velocity field for the grains and 7? c d ~ 1, chondrules 
will sweep up the dust before significant aggregation takes 
place and adust ~ ao; contrarily, if Brownian motion domi- 
nates (or when r cn » f s ), the monomers will first collide with 
each other before being accreted by chondrules. 

The question that remains is what this means for the poros- 
ity of the rim. Assuming a hit-and-stick process, in which the 
accreting dust particles are all of equal size and much smaller 
than the chondrule/compound, the structure of the rim will re- 
semble that of particle-cluster agglomeration (PCA). Thus, if 
the dust particles are solid monomers the filling factor of the 
rim indeed equals (ppcA- On the other hand, if the accreting 
dust particles are fluffy aggregates, but still smaller than the 
chondrule, the resulting filling factor of the rim will be less 
than </>pcA- The precise filling factor will be determined by the 
amount of interlocking between the aggregates but as a crude 
upper limit the aggregated may be approximated as a homo- 
geneous porous sphere such that the packing process of the 
dust aggregates is PCA. Then, the filling factor of the porous 
dust, p d, is just the product of 0pca (caused by hit-and-stick 
packing) and the porosity the dust aggregates already contain 

(0dust), i-e., <f>pd = ^pcA^dust- 

In yet another collisional growth scenario we envision that 
chondrules are mixed into a dust region after this dust has 
aggrega ted into dust balls but before the dust balls are com- 
pacted (lOrmel et al.l 120071) . In any case we assume here 



that the dust consists of porous aggregates and parameter- 
ize its filling factor by the pd = ^pcA^dust relation. Using 
equation d2Tb , ad us t is the parameter that regulates the fluffi- 
ness of the dust accretion process and we run models at dif- 
ferent values of Adust t° tes t its importance and sensitivity. 

2.4.3. Compound-compound accretion 

The timescale for chondrule-chondrule accretion, f cc , is 
similar to the dust-chondrule timescale, equation dZTl i. Dur- 
ing the collision, a fraction of the porous dust is lost to the 
compact dust phase. Since the porous phase filling factor is 
always less than that of the compact phase there is always a 
net loss in geometrical volume when two compounds collides, 
i.e., V < V\ + V-2. This, we call the 'conservative approach' 
(no fractal accretion of compounds). Alternatively, collisions 
of compounds (consisting of one or more dust-rimmed chon- 
drules) may be in the hit-and-stick regime. This would occur 
if the impact energy is absorbed locally and is not communi- 
cated to other parts of the compo und. In that case t he com- 
pound packing proceeds fractally. Orme l et ail d2007h provide 
an expression for the growth of V in the hit-and-stick case for 
particles of different size, derived by an interpolation from the 
PCA and CCA limiting cases, i.e., 



-3/4 



V = V, 



(28) 



where V\ is the volume of the largest of the two particles that 
meet and 5 ^ 0.95. The growth of the porous phase V pd 
then results from the gain in V through equation d20l >. The 
porous phase is then a mixture of porous dust and voids and 
the geometrical volume becomes a balance between hit-and- 
stick packing of chondrules (increasing V) and compaction of 
porous dust (decreasing V). This, contrary to the conserva- 
tive approach in which the total volume always decreases at 
collision. 

For example, if V2 = j^Vi, equation d28| i shows a volume 
0.045 Vi is added to the porous phase, decreasing its filling 
factor. It then depends on / comp how the net growth of the 
porous phase turns out. In the initial stages of coagulation 
/comp is often very low and, therefore, fractal accretion of 
compounds can become very important in enhancing the geo- 
metrical volume of the compounds. 

Although fractal accretion of compounds increases the vol- 
ume of the porous phase, no mass is transferred to it. It is 
only the filling factor that is affected, in its turn affecting the 
aerodynamic properties of the compound. Eventually, due to 
compaction, all models run out of porous dust and the dust in- 
side the final objects - whether fractal accretion is involved or 
not - has the same filling factor c d = 0.33. The lower filling 
factor of the porous phase during the collision process merely 
reflects the voids between the chondrules that are created in 
the models with fractal accretion. As velocities increase, how- 
ever, the fractal structure must collapse. We suspect, further- 
more, that structures of very low filling factor are too weak to 
surviv e the more violent collisions (see, e.g., iParaskov et alj 
l2007h . 

2.5. Collisional scenario 

We have proposed a model where chondrules - in the pres- 
ence of dust particles - acquire rims of fine-grained dust, 
which help them stick together, and discussed various colli- 
sional scenarios for this growth process. Here we briefly sum- 
marize the envisioned scenario from a chronological view- 
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point, and emphasize the (free) parameters the model con- 
tains. 

At the start of our simulation (f = 0) a population of chon- 
drules encounters a reservoir of dust particles of characteris- 
tic size adust (see {32.4.21) . The physical conditions of the disk 
(e.g., gas density px, turbulent strength a; see {32.2b determine 
the relative velocities between the chondrules/compounds. 
Another important parameter is the spatial density of dust (pd) 
or, rather, the dust to chondrule density ratio, 7? c d, since this 
determines the thickness of the rims. These, and other phys- 
ical conditions at the start of the simulation determine the 
subsequent accretion process. First, chondrules start to ac- 
crete the dust aggregates (see {32. 4. 21 >. resulting in a porous 
rim of filling factor, p d. When rimmed chondrules collide, 
part of this porous structure collapses to c d = 0.33 filling 
factor through the initiation of rolling motions. This dis- 
sipates the collisional energy and, if enough porous dust is 
present by the criterion of equation (fTol l. the two chondrules 
stick and a compound is created. In this way compounds are 
created and many chondrules can be stuck together until the 
accretion process is terminated when both the amount of free- 
floating dust and the porous dust inside the compounds have 
become insignificant. The end product is an inert population 
of compact-dust rimmed chondrules and compounds that only 
bounces. Collisional fragmentation is not explicitly included 
in the model, but we can a posteriori compare the velocities 
with a critical threshold (~ m s _1 ) to verify its importance. 

3. MONTE CARLO COAGULATION 

The physical model of chondrule accretion contains many 
free (i.e., unknown) parameters. In a statistical study of com- 
pound coagulation, we will sample these free parameters at 
discrete intervals such that a grid of models is created (see 
In this section we briefly describe the Monte Carlo coag- 
ulation code used to calculate the collisional evolution of the 
compounds. 

Compounds are characterized by four numbers - e.g., m, 
/d and / p to determine the mass inside each phase, and p d 
for the filling factor of the porous phase. Therefore, a Monte 
Carlo code, rather than the multi-variate Smoluchowski equa- 
tion, is the obvious way to solve the collisional evolution. In 
our code we do not keep track of the individual positions of 
eac h constituent unit (th e monomers) within a compound as in 
e.g jKempf et al.l(l 19991) but identify each compound by these 
four numbers. In this way we have a good characterization of 
the internal structure of the compound, although the precise 
internal structure cannot be retrieved. 

The code we use is called event driven; i.e., the timestep 
At de f ines the time between two con secutive events (Gillespie 
119751: ISpougelll985t IUffmanlll991l) . Here, events are colli- 
sions between two compounds (see below for the dust). The 
collision rate C,y gives the probability of a collision - CijAt is 
the probability of collision in the next (infinitesimal) timestep 
At involving compounds i and j, i.e., 

<t Au N N 

Ci >-^r 1: and c t°t = 2Z c '> (29) 

i j>i 

in which cry is the collision cross section 5 between particles i 
and j, Avij the relative velocity between the two compounds, 

5 In general one must distinguish between the collisional cross section 
(which gives the reaction rate of the two species) and the geometrical cross 
section (which is the average projected area of the particle that determines the 
coupling to the gas). Here we will sim ply equate them as in erg = n(ai + a,) 
and ignore the small discrepancy (see lKrause & Blum 2004). 



and N the total number of particles in the simulation. The 
volume of the simulation, r V, is determined from the spatial 
density in chondrules, the constant p c , and the total mass in 
chondrules, i.e., *V = Yli=\ _ /d,;)/Pc- From these quanti- 
ties the timestep is defined by At = -C t ~| In 7, with r a random 
deviate. The particles that are involved in the collisions are 
also determined randomly, weighted by their collision rates 
Cjj. Then, using the recipes outlined in ^2] the outcome of 
the collision - sticking or bouncing - is determined. In either 
case, the parameters of the new or modified compounds are 
re-computed. (In the case of bouncing the change is reflected 
in a smaller size, a, due to the compaction.) Subsequently, the 
new collision rates of the particles (i.e., the {Ca} and, if the 
second particle due to bouncing is still present, the \Cjk) for 
k — 1 ... N and k + i, j) are re-computed. These updates of 
the collision rates are the most CPU-intensive part of the code. 
With it one cycle is completed, after which a new stepsize At 
is determined and the steps repeat themselves. 

Equation d29i i involves the total relative velocity between 
the particles. To calculate Auy we use thermal, turbulent and 
systematic velocities (eqs. J4), (8) and ifTOl ). adding them up 
in quadrature. Strictly speaking, the zero-dimensional nature 
of the MC-model is inconsistent with dispersal of particles 
(particles do not have a positions); however, when the drift 
is modest the change in the physical environment is negli- 
gible and we can still use the MC-approach. For the radial 
drift this assumption applies only when the total drift is small 
compared to the initial location of the particle, i.e., AR <c R, 
such that the same physical conditions apply throughout the 
simulation. We will a posteriori check whether radial drift is 
significant. 

Apart from collisions between compounds, we also keep 
track of dust accretion. This is, however, not implemented in 
a Monte Carlo fashion: it would have made the code very slow 
since the tiny dust particles far outnumber the chondrules. In- 
stead, we catalogue the cumulative dust mass that is accreted 
by the compounds over the timesteps; i.e., for compound i we 
increase the amount by najAu^pdi^At. Only when this mass 
exceeds a certain fraction (say / llp d = 10~ 3 ) of the total mass 
of the compound, this quantity is added as porous dust to the 
compound and the /d,/ p parameters as well as the collision 
rates are updated (eq. |29l ). Although this procedure makes 
the program still a bit slow at the initial stage of the simula- 
tion, it is definitely much faster than updating all ^N(N - 1) 
collision rates at every timestep. The (decreasing) amount of 
free-floating dust, pd(t), is computed in this way. We have 
examined the sensitivity of this mechanism on / up d and found 
that / up d = 10~ 3 is accurate, while much more efficient (faster) 
than, e.g., 10~ 5 . 

The strong point of the MC code is that it can deal with 
many structural parameters and that it is transparent and 
straightforward; the weak point, however, is its low numer- 
ical resolution. Given the complexity of the model and the 
large number of models we intend to run, the number of par- 
ticles (AO we use in the simulations is a few thousands at 
most. To prevent the resolution from deteriorating (a colli- 
sion resulting in sticking decreases N by one) we artificially 
stabilize the total particle number by a procedure called dupli- 
cation. In this process, one particle is r andomly chosen and 
dupli cated from the existing population ( Smith & Matsoukas 
1998). Subsequently, *V is increased proportionately such 
that the total density in chondrules, p c , stays constant. This 
procedure is called the constant-A algorithm - an algorithm 
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far su perior in terms of accurac y to the constant "V algo- 
rithm ( Smith & Matsoukai 1 19981) . and we have previously 
shown that it is able to calculate large ord ers of growth, espe - 
cially when the size distribution is narrow (lOrmel et al.l2007l) . 
Through the duplication mechanism, furthermore, a distinc- 
tion can be made between 'duplicates' and 'distinct species,' 
and it is actually the latter that we keep constant, such that the 
total number of compounds involved can be much larger than 
a few thousands, also improving the efficiency of the model. 

4. RESULTS 

In our models we generally recognize three stages in the 
growth process: hit-and-stick dust accretion (increasing the 
porosity), compound accretion (i.e., growth), and compaction 
with accompanied stalling of the growth. The balance be- 
tween these phases controls the size of the resulting com- 
pounds, while their relative importance and 'timing' are deter- 
mined by the adopted model parameters. In §14. II we discuss 
two illustrative cases, focusing on the temporal stages during 
their evolution. Then, in {34.2l we investigate the sensitivity of 
the other parameters by means of a parameter study. 

4.1. Individual model runs 

Figures [6] and [7] show detailed results for two individual 
runs of the simulation with (default) parameters of gas den- 
sity px = 1, Adust = 10 /mi, 7? gc = 100 and = 1, 
y=19 erg cm' 2 at a distance of 3 AU (see Tabled- In these 
figures panels A correspond to a model with a = 10~ 4 , while 
a = 10~ 6 in panels B. In Fig.|6]several (mass-averaged) quan- 
tities are shown as function of time, while in Fig. [7] the size 
distributions of compounds are shown at three points during 
their evolution: i) t — (the initial size distribution of chon- 
drules); ii) the time at which 50% of the dust is accreted; and 
Hi) the time at which a negligible amount of porous dust re- 
mains (the final distribution). The negligible criterion is met 
when both the porous dust mass within all compounds as well 
as the density of free-floating dust are less than 0.1% of the 
initial dust mass. 

In Fig. [6] we make use of mass-weighted averages. For ex- 
ample, the mass-weighted average size of the population is 
defined as 

(a) m = — , (30) 

£; m 

where the summation is over all particles of the simulation. It 
gives the mean size in which most of the mass of the popu- 
lation resides, and is more appropriate to describe the popu- 
lation than the average size, (a). In particular, adding a large 
number of small particles with negligible mass (density) to 
the population, decreases (a) but leaves (a) m unaffected. In 
the following the prefix 'mw-' is used as an abbreviation for 
'mass-weighted average of the distribution.' 

Figure|6j\ shows that dust is accreted on times cales of a few 
10 2 yr, which agrees well with previous studies (Cuzzi 2004). 
This causes the size distribution (Fig. [7]) to shift towards larger 
sizes: the accretion of porous dust particles at low filling fac- 
tor significantly increases the geometrical size of the com- 
pounds. The dotted curve in Fig.[6}\ shows the mw-filling 
factor of the accreted dust. At the start of the simulation 
this equals p d = 0dust0PCA; however, collisions are energetic 
enough to compact the porous dust on a global scale. The de- 
crease in friction time (solid grey curve), caused by the accre- 
tion of porous dust, therefore is only modest. (Note that even 
accretion of <\> — c d = 33% dust on chondrules would cause 



the friction time to decrease). Consequently, the sticking rate 
never increases much beyond ~ 50%. After ~ 10 2 yr the ac- 
cretion of porous dust cannot keep pace with the compaction 
and sticking fails, resulting in a decrease of the compounds 
geometrical size (dashed curve). This results in a 'retrograde 
motion' of the final size distribution curve in Fig. [7j\. In the 
a = 10~ 6 model (Fig. [6}}) the collision velocities are much 
lower and, different from the a = 10 4 model, the porous dust 
does not experience compaction for a long time. The sticking 
rate then increases to almost 100%. However, depletion of 
dust triggers the end of the growth phase; growth is quickly 
terminated by the mutually enforcing processes of rim com- 
paction and increasing velocities. From these panels it is clear 
that much growth can be achieved when relative velocities are 
kept low during the dust and chondrule accretion. 

Although the Monte Carlo code does not keep track of the 
position or size distribution of chondrules within compounds, 
we can still extract useful statistical information from the 
model runs. One such statistic is the distribution of the dust 
over compounds of different size: is the dust primarily ac- 
creted by individual chondrules or by large compounds con- 
taining many chondrules? The results are presented in the 
histogram of Fig. [8] The x-axis denotes the number of chon- 
drules a compound contains (N) and is divided into logarith- 
mic bins of base 2, i.e., the first bin corresponds to single- 
chondrule compounds, the second to compounds that contain 
2 or 3 chondrules, the third to 4 - 7 chondrules, etc. The 
i/-axis gives the relative fraction of the dust that first ac- 
cretes onto a chondrule or compound with size in each bin; 
that is, Fig.[8]only reflects the dust accretion history and does 
not include the subsequent re-distribution of dust due to co- 
alescence of compounds (which would shift the dust-rimmed 
chondrules to a larger compound bin). The relatively high 
level of the first bin (single-chondrule compounds) reflects 
dust accretion during the early phase of the simulation where 
individual chondrules provide a high surface area and the den- 
sity of free-floating dust is highest. In simulations with strong 
turbulence this fraction becomes very high: the dust is then 
accreted by single chondrules only. But even in the case of 
low a single chondrules are responsible for a significant share 
of the dust sweep-up, as the a = 10~ 6 results show. Besides, 
larger compounds also have a larger surface to spread this dust 
over; rims created by dust accretion are therefore thickest on 
chondrules. 

4.2. Parameter study 

Figure|9] presents the results of the parameter study. The 
free parameters (Table [5]) are distributed over a grid such that 
each grid point corresponds to a unique model. In total a few 
thousand distinct models are run. Each model is run a few 
times to account for stochastic effects in their results (typi- 
cally ~10% or so). For each free parameter the models are 
ordered by the grid-values of the parameter, corresponding to 
the panels in Fig. [9] Two output values are shown: the ra- 
tio of final mw-radius to initial mw-radius, or growth factor 
(crosses) and the time at which the dust is depleted and the 
simulation terminated, or simulation time (diamonds) (they 
share the same i/-axis; for clarity the timescale error-bars are 
slightly offset in the ^-direction). The symbols denote the log- 
arithmic averages of all models at the grid-values and the error 
bars indicate the range in which 50% of the models fall. This 
spread can be huge since it is primarily determined by the 
spread in the other parameters (and therefore nowhere close 
to Gaussian). The same holds for the averages: these can be 
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Fig. 6. — (A) A run of the compound accretion model with a — 10 , px — 1. ^dust = 10 jum, y = 19 erg cm , K = 3 AU, 7? gc = 100 and "R c( j = 1 (the default 
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with a- = IO" 6 . 
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Fig. 7. — Size distributions of compounds corresponding to the runs in Fig.|6]for the a = IO" 4 model (A) and the a = 10-° model (B). Shown are the initial 
distribution (crosses), the distribution at the time where 50% of the dust has been accreted (plus-signs) and the final distribution (diamonds). Note that compaction 
has the effect of shifting the distribution to the left. 

prisingly irrelevant. For example, the dependence on the size 
of the dust fluff-balls, adust, defining their porosity (eq. J2T|), 
is only modest (panel B), and also the latter two panels do 
not show clear trends. Panels H and I directly give the ratio 
between the two data points, instead of the exponent of the 
power-law fit. Panel H shows the effects of taking the geom- 
etry of the collision into account. / geo = X where < X < 1 
means that / geo i s dete rmined by the geometry of the collision 
as discussed in §12.3.31 whereas / geo = 1 indicates that all dust 
is available for compaction. However, allowing for a lower 
/geo also reduces the maximum amount of dust that is com- 
pacted. Apparently, these two effects largely compensate. A 
similar insensitiveness is shown in panel I: whether we allow 
for fractal accretion of compounds ('1') or not ('0') does not, 



arbitrarily scaled up or down by giving more weight to ex- 
treme models in the parameter study. 

However, the value of Fig. [9] lies not in its absolute num- 
bers but in the trends that emerge from the parameter varia- 
tion. The lines indicate this trend and their slopes are given 
in each panel. From these, it is seen that timescales are pri- 
marily determined by turbulent a (velocities), nebula location 
R (densities) and the chondrule density (panel E). Parame- 
ters that favor large growth of compounds are low a (panel 
A), high gas densities (panel C), low chondrule-to-dust ra- 
tios (panel F), and high surface energy densities (panel G). 
Growth is favored in these models due to the moderate rel- 
ative velocities (panels A and C) or better sticking capabili- 
ties (panels F and G). Other parameters are sometimes sur- 
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Fig. 8. — Relative fraction of dust accreted by compounds of different size 
for the a = 10~ 4 and a = 1CT 6 simulations. The compounds are placed in 
bins according to the number of chondrules they contain. The bins are expo- 
nentially distributed by factors of two. The histogram shows the distribution 
of the dust over the compound sizes (in terms of number of chondrules inside 
the compound) at the time of the dust accretion. Single chondrules (first bin) 
accrete a significant fraction of the dust. 



TABLE 5 

FREE MODEL PARAMETERS 



(1) 


(2) 


(3) 


(4) 


(5) 


turbulent strength 


a 




5 


per 7 — i(T 3 ] 


size of dust aggregates 


^dust 


cm 


3 


[lO^ 4 — 10" 2 ] 


gas density" 


Px 


g cm -3 


3 


[1—100] 


nebula location 


R 


AU 


3 


[1, 3, 10] 


gas-chondrule ratio 


Kg C 




2 


[10, 100] 


chondrule-dust ratio 






2 


[1, 10] 


dust composition b 


y 


ergs cirT 2 


2 


[19, 370] 


compaction mode 


/geo 




2 


[X, 1] 


fractal accretion of compounds 




2 


[yes, no] 



Note. — List of free model parameters. Columns denote: (1) pa- 
rameter description; (2) symbol; (3) unit; (4) number of grid points per 
parameter; (5) parameter range, with a grid point at every factor of 10, 
unless otherwise indicated. See also Table |2]for other (fixed) parame- 
ters. 

a px'- gas density over MSN model at 3 AU. 

b y: energy surface density. The high y model corresponds to ice as the 
sticking agent (10 AU models only). 

c /geo = 'X': /geo is computed after the procedure outlined in 32.3.31 
/geo = 1 '• use /geo = 1 always. 

d Whether a hit-a nd-stic k packing model for compounds (leading to a 
fractal structure, 32.4.3) is adopted or not. 

on average, make a difference. Note, however, that the values 
shown in the panels are averages; in some individual simula- 
tions we do see a notable increase when fractal accretion is 
turned on. 

Panel A shows that the positive correlation between growth 
factor and decreasing a breaks down for values below a < 
10~ 5 : the growth flattens out and reaches a constant level. The 
reason is that for low values of a, and, subsequently, large 
compounds, radial drift motions quickly take over from tur- 
bulent motions such that the evolution becomes insensitive to 
a. The high radial drift obtained when particles approach the 
St = 1 regime forms a barrier for further growth. Low a com- 
bined with high gas densities delay this transition since com- 
pounds are now much better coupled to the gas (lower Stokes 



numbers) meaning much growth early on. Yet, Stokes num- 
bers inevitably grow to values near unity, and in most cases 
the resulting rjv^ drift velocities (eq. (8l) stall growth below 1 
meter. 

Panel D shows that growth depends only modestly on neb- 
ula radius, R. Here, the positive correlation with growth factor 
is a bias resulting from the high y 'ice models' (panel G) - ice 
promotes sticking - that are only present at R = 10 AU. Thus, 
despite the fact that almost all nebula parameters scale with 
R, their combined effect does not result in a clear trend that 
favors growth. For example, larger nebula radii mean lower 
densities and higher Stokes number, increasing the velocity in 
the initial stages, but this is offset by a (slightly) lower sound 
speed, and the better sticking agents that are available. 

In Fig. [10] all models are combined in a scatter plot of to- 
tal (mass-weighted) radial drift against the final mw-size ob- 
tained in the simulation. The few models that cluster around 
the meter size are all ice models (y = 370 and R — 10 AU, 
indicated by triangles). Some of them do manage to cross 
the (St)„, = 1 barrier (at 10 AU and px — 1 this already 
happens at a few centimeters) but do not make the jump to 
planetesimal sizes. In models that during their growth drift 
less than AR ~ R the local assumption is justified; for models 
that drift over several AU-distances, however, the approxima- 
tion we used in the calculation of the collisional evolution, 
i.e., that the physical conditions stay the same, breaks down. 
Note, however, that the drift distances in Fig. [10] are upper 
limits: radial drift slows down with decreasing R due to a 
better coupling to the gas, or may diminish when collective 
effects become important (see 35. It . 

4.3. Importance of fragmentation 

The models neglect the possibility that high velocity colli- 
sions will fragme nt, rath er than merely compact or bounce, 
compounds. In 32.3.51 it was estimated that at ~ m s _1 
velocities, fragmentation becomes likely, starting with ero- 
sion, followed by catastrophic disruption of the compound. 
This threshold can now be compared to the maximum veloc- 
ities attained at a particular Stokes number (see Fig. [TJ, i.e., 
~ /7^ 2 c g St and ~ c'^CgSt 1 ^ 2 for the systematic drift and tur- 
bulent velocities, respectively. Inserting the final mw-Stokes 
number into these expressions, we obtain a criterion whether 
fragmentation is of importance. As the critical velocities we 
take 2 m s _1 for turbulence and 6 m s _I for radial drift (see 
Fig- El- It can be shown that for systematic velocities scaling 
proportional to size the collisional energy peaks at size ratios 
of e ~ 0.5, which, according to Fig. [4] corresponds to a frag- 
mentation velocity of ~ 3 m s _1 , or a ~ 6 m s _1 radial drift 
velocity for the largest particle. For ice models (triangles) 
the fragmentation threshold is increased by another factor of 
four, reflecting their higher y. In Fig.[10]the models that have 
crossed the threshold velocity are indicated by an open square, 
whereas black squares indicate velocities that stay below the 
threshold. The grey squares are models in which collective 
effects could have had a significant reduction in relative ve- 
locities and drift rates, due to concentration of compounds 
near the midplane (see below, 35.11 ). However, these sub- 
tleties are presently not taken into account in the simulation 
and it remains unclear whether fragmentation is an important 
phenomenon in models for which settling is important. 

For most models in Fig.fTOl fragmentation is not a seri- 
ous concern. This is a natural result as compaction precedes 
fragmentation and growth stalls before reaching the fragmen- 
tation threshold. However, for large, fluffy compounds the 
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a drift of 1 AU. The dashed horizontal lines indicate compound sizes of 1 cm and 1 m, respectively. 



compaction is more pronounced, resulting in a significant de- 
crease in surface area-to-mass ratio, increasing the Stokes, 
and thereby possibly breaching the threshold for fragmenta- 
tion. Also note the fragmentation models (open symbols) at 
the bottom of Fig. [10] in these the fragmentation threshold 
was already exceeded at the start of the simulation. 

Figure QT| shows a small subset of models from Fig. [TOl that 
takes away the redundancy (caused by less influential param- 
eters) and focuses on the more plausible scenarios. More 
specifically, Fig. QT| shows models limited to the following 
parameters: ??g C = 100, Red = 1.0, adust = 10~ 3 cm, px = 
1 or 10, R — 3 AU, y — 19 erg cirT 2 ; furthermore, we assume 
only local compaction (/ geo = X) and assume collisions be- 
tween compounds are not in the hit-and-stick regime ( jj2.4.31 l. 
Only 12 models are then shown with the only free parame- 
ters being a (all 6 distinct values) and px (2 values). Table 
[6] shows various output values corresponding to the 'top ten' 
models of Fig.Qj] ordered after final mw-size; for example, 



the maximum velocities due to systematic and turbulent mo- 
tions. This shows that for these low-a models systematic drift 
velocities (Col. (7)) quickly become dominant over turbulent 
motions (Col. (8)). 

5. DISCUSSION 
5.1. Collective effects in a settled layer 

Despite the ability of the chondrule-sticking model to tweak 
many parameters to optimize the growth, compounds never 
achieve planetesimal sizes. Ultimately, m/s or higher veloci- 
ties are unavoidable in any model due to the radial drift; that 
is, compounds inevitably reach (and have to cross) the St = 1 
bottleneck at which their radial drift velocities peak. The stud- 
ied accretion mechanism - chondrule sticking by compaction 
of initially fluffy dust - is simply too weak to grow past the 
St = 1 bottleneck. 

There is one issue, however, that is unaccounted for within 
the framework in which the simulations are performed: if the 
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TABLE 6 

DETAtLED RESULTS 



a 


PX 


(a) m 
[cm] 


<st> m 




Av tmb 

[cm s _1 ] 


A „sys 
[cm s _1 ] 


A[ ,sy S ,CE 

[cm s _1 ] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


io- 7 


10 


8.5 


2.6 x 10~ 2 


4.9 


5.9 


358.6 


10.3 


10" 6 


10 


8.0 


2.5 x 10~ 2 


1.5 


18.0 


336.9 


54.6 


IO" 5 


10 


4.8 


1.5 x 10~ 2 


0.4 


44.3 


204.4 


108.2 


10- 7 


1 


0.8 


2.5 x 10~ 2 


4.8 


5.8 


345.0 


10.2 


IO" 6 


1 


0.7 


2.3 X 10~ 2 


1.4 


17.4 


316.7 


53.4 


10" 4 


10 


0.6 


2.0 x 10~ 3 


0.0 


51.2 


27.4 


25.1 


IO" 5 


1 


0.5 


1.4 x 10~ 2 


0.4 


42.8 


190.9 


103.8 


io- 3 


10 


0.1 


3.9 x IO" 4 


0.0 


71.4 


5.3 


5.2 


10~ 4 


1 


0.1 


4.0 x IO -3 


0.1 


71.9 


53.9 


47.9 


IO" 3 


1 


0.1 


3.3 x IO" 3 


0.0 


206.7 


44.5 


42.8 


io- 2 


1 


0.1 


3.2 x 10~ 3 


0.0 


641.7 


42.9 


42.0 


io- 2 


10 


0.1 


3.5 x 10~ 4 


0.0 


212.3 


4.7 


4.6 



Note. — Detailed results from 12 selected models (see text), ordered after 
final mw-size, Col. (3). The columns denote: (1) turbulent-a; (2) gas-density 
enhancement (restricted to 1 or 10); (3) final mw-size; (4) final mw-Stokes 
number; (5) final midplane dust-gas density ratio would settling have been in- 
cluded; (6) turbulent velocity contribution, or'' 2 St'' 2 Cg; (7) systematic drift after 
equation (8); (8) systematic drift due to collective effects after equation 13 It with 
Col. (4) for St. 

turbulence is weak enough, in addition to moving radially, 
compounds can also settle into a dense layer at the midplane 
as their Stokes numbers increase. When the density of solids 
at the midplane exceeds the gas density, the gas is dragged 
with the particles (instead of the other way around), result- 
ing in gas velocities that tend to become closer to Keplerian, 
which subsequently diminishes the radial dri ft and relative ve- 
locities of particles. Nakaga wa et alJ 0986) have solved the 
equations of motion in such a two-fluid medium analytically 
for a single particle size (or Stokes number); the radial drift 
velocity now becomes (instead of eq. (8)) 



2St 



(31) 



St 2 + (1 +p p (z)/p g ) 2 ' 
where p p (z) is the total density of particles at a height z 



above the midplane. For a generalized solutio n over a particle 
size di stribution see Weidenschilling (1997) or Tanak aet al.l 
(2005). Thus, in a dust-dominated layer the radial drift of in- 
dividual particles depends through p p (z) on the density of all 
other particles: a collective effect. The particle concentration 
can be found by balancing the gravitationally induced settling 
rate with the diffusion rate of a particle, assuming a steady- 
state distribution. The scaleheight of the resulting part icle 
distribution, h p , can be calculated as (Dubr uTle et al.ll 19951) 

h v = -=J=, (32) 

^^TTs 

where S = St/ a. Under conditions of initial cosmic abun- 
dances, in order to reach p p ~ p g the particles must settle 
into a layer of thickn ess one-hundredth of the gas scaleheight, 
requiring S > 10 4 (ICuzzi et al.1 120051) . This may occur for 
chondrules in very low-cr environments, or, at moderate a, 
only for large compounds during their growth and settling 
stage. In Table|6]we have calculated the density enhancement 
(Pp/Pg, Col. (5)) and the corresponding velocities (Au sys ' CE , 
Col. (8)) for a few selected models at the end of their simu- 
lation, where we fixed most parameters at their default 3 AU 
values, except for a and px- Note that in the context of our 
current model setup collective effects are purely hypothetical 
(we treat p 9 /p c = fi g c as a constant); the columns of Table[6] 
therefore merely provide an indication of what could be ex- 
pected had settling-effects been included. In these calcula- 
tions we have used the mass-averaged Stokes number of the 
population (Col. (4) of Tabled as the Stokes number that en- 
ters equations PTT ) and d32b . The last two columns of Table [6] 
show that collective effects (p P /p g > 1) quickly reduce the ra- 
dial drift. In a future study, we intend to investigate the effects 
of the particle concentration on the compounds' growth. 

There is yet another subtlety involved when collective ef- 
fects (i.e., a dust-dominated midplane) be come important. 
This is the Kelvin-Helmholtz instability ( We idenschilling! 
1980), caused by the shear between the two fluids now mov- 
ing at a relative velocity of AV, the azimuthal velocity differ- 
ence between the gas in the particle-dominated and the gas- 
dominated layer. For shear turbulence the turbulent viscos- 
ity is v T ~ (AV) 2 /QRe* 2 jCuzzi et alJI 1993b . where Re* is a 
critical Reynold s number at which the flow starts to become 
turbulent, which lDobrovol skis et alj (1 19991) found to be Re* ~ 
20 - 30. Also, the large eddy turnover frequency in shear 
turbulence (Q e ) can become much larg er than Q, depending 
on the thickness of the shear layer (see Weidenschillina l2006l 
for how Q e depends on the particle density structure, p p (z)). 
Equating v T ~ (AV) 2 /QRe* 2 with (u L hear ) 2 /Q e then provides 
the expression for the shear turbulent (large eddy) velocity, 



..shear 



, shear 



1/2 



0.033 (52 



1/2 



AV < 0.18/7 1/2 - 



Co, 



Cl) Re* 

(33) 

where the upper limit assumes (Q e /£T) ~ Re* = 30 and 
AV = rjvk = rj l t 2 c g . This corresponds to the situation where 
the shear layer is thin (meter-size or larger particles; the 
shear layer cannot become thinner than the Eckman layer, 
see ICuzzi et"ai1ll993l) . In that case, setting y* ear = a'^ ar c g 
the equivalent a value for shear turbulence becomes a s h e ar ~ 
3 x IO -5 . This is an upper limit; for smaller particles, or a 
size-distribution of particles, both Q e and AV are lower and 
tfshear decreases as well. Shear turbulence may therefore be 
much more conducive to compound growth. 



Accretion of chondrules in the solar nebula 



19 



Future studies must show whethe r these effects enable 
growth to planetesimal sizes. Recently, Johansen et al. (2007) 
have suggested that concentration of meter-size particles (St ~ 
1) in certain azimuthally-oriented near-mi dplane high pres- 
sure zones, which form between large turbulent eddies, might 
lead to gravitationally bound clumps with the mass of plan- 
etesimal size objects. The results from Johansen et alj (120071) 
were most pronounced when the turbulent intensities were 
moderately high (this leads to the largest radial pressure con- 
trast), suggesting values of a ~ 10~ 3 . However, our results 
suggest that it is difficult to grow a population of meter-size 
boulders in the first place under such conditions. The maxi- 
mum growth (in terms of Stokes number) our models achieve 
for a = 10~ 3 is St ~ 5 x 10~ 3 at 3 AU (essentially no growth 
at all: just dust-rimmed chondrules). 10 AU ice models do 
somewhat better: St ~ 7.4 x 10~ 2 . Even if they can form, a 
population of meter-sized boulders may be difficult to main- 
tain if these originated from dust-coated, solid chondrules 
as modeled in this paper. In §12.3.51 we estimated that frag- 
mentation occurred at a critical velocity of ~ 10 2 cm s _1 , 30 
times smaller than the expected value of St = 1 particles for 
a = 10~ 3 . This translates into a specific kinetic energy for 
disrupt ion of Q* = 10 4 erg g 1 , much lower than the criti- 
cal Q* Johanse n et al.l (120071) adopt (for aggregates of solid 
basalt objects, as taken from iBeria 12.0001) . Thus, our results 
indicate that it may b e difficult for the instability described by 
iJohansen et al.l (120071) to become viable in the turbulent inner 
(ice-free) nebula. 

In the outer solar system, however, conditions may be more 
favorable to growth in a turbulent environment. First, if 
ice acts as the sticking agent Q* may be over an order of 
magnitude larger, reflecting the scaling with the surface en- 
ergy density parameter, y. Second, if chondrule formation is 
not common the particles grow directly from aggregates of 
tiny grains to larger aggregates and therefore contain roughly 
twice as much mass in small grains as dust-rimmed solid 
chondrules. Moreover, Stokes numbers for the same parti- 
cles increase with larger heliocentric radii (R) due to the lower 
gas densities. Therefore, at large R the St ~ 1 regime is 
reached at smaller sizes (centimeters), which may be some- 
what more difficult to disrupt (i.e . , higher Q *) than m-size 
bodies dHousen & Holsapplelll990T:lBenzll2000l) . (In our sim- 
ple estimate of Q* we do not have a size dependence, though.) 
Still, it is hard to see that > 10 m s _1 velocity collisions be- 
tween equally-sized particles, even under these most favor- 
able conditions, will not result in disruption; but this should 
of course really be tested by experiments. 

On the other hand, it also seems sensible to pursue incre- 
mental growth scenarios which take place in quiescent (or 
low-a) nebulae. Due to the relatively low effective a-values 
for shear turbulence derived above, it is the radial drift mo- 
tions that will provide the limits to growth. However, even a 
modest reduction of radial drift motion by a few factors due 
to collective effects may already be sufficient to pr event catas- 
troph i c collisions as particles reach St — 1 (see, ICuzzi et all 
1993; Weidenschilling 2006). Recall that fragmentation is 
easiest for nearly equal-size particles, which collide at very 
low velocities due to their systematic, nonturbulent motions. 
Moreover, the compound size distribution, which is a function 
of height, also determines whether collisions are beneficial 
to the growth; for example, if the size distribution is nearly 
monodisperse this will certainly favor growth in the nontur- 
bulent cases. Since all these effects will vary with height, 
however, it is difficult to predict how these effects will unfold, 
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Fig. 12. — A sketch of a cross-cut through a compound at the final state of 
our model. The compound contains two phases, present in equal proportion 
by mass: chondrules (black) and <j> = 0.33 compact dust (grey). The cross- 
cut introduces a selection effect and shifts the chondrule size distribution to 
bigger chondrules. The chondrules are placed at random but a certain distance 
between the chondrules is preserved (see text). 

and which parameters are key. Clearly, additional modeling 
is needed, where we may even combine these two different 
modes of turbulence since it is quite natural to expect t hat dif- 
feren t physical processes operate at different heights (ICieslal 
120071) . However, incremental growth in the dense, particle- 
dominated midpla nes of nonturbulent models then proceeds 
extremely rapidly dCuzzi et alj|1993t IWeidenschiHing; l2000l) 
which is contrary to the ev id ence from meteorites and aster- 
oids fsee lCuzzi et~ai1l2005l or lCuzzi & Weidenschillingll2006l 
for a discussion). 

5.2. Dust rim and matrix 

Figure Q~2] provides an illustration of the internal structure 
of the objects obtained at the end of our simulation. In Fig. Q~2] 
it is assumed that each chondrule (black circles) is surrounded 
by a dust rim at least ~ 40% of the chondrule's mass, corre- 
sponding to the amount of dust accreted by individual chon- 
drules (see Fig. [8j. This translates into an outer rim radius 
that is a factor of 1 .3 larger than that of the chondrule and is 
reflected in the inter-chondrule spacing of Fig. [12] For the 
remainder the chondrules are positioned at random. (Note 
that Fig. [12] does not follow directly from the Monte Carlo 
collision model since we do not keep track of the positions 
of chondrules within compounds and cannot 'reconstruct' a 
compound.) Furthermore, we have assumed the initial chon- 
drule size distribution holds, but accounted for the bias a cros s 
cut introduces to the observed structure (Eisenhour 1996). 
The chondrule size distribution is therefore skewed towards 
larger chondrules. 

However, the picture of a uniform 67% porosity dust phase 
between chondrules contrasts with the meteoritic record. 
Here, the dust is compacted to a much larger extent and - 
at least in the pristine CM chondrites but in others as well 
- can be divided into two dis tinct components: the fine- 
grained, low porosity (10-20%; Trigo-Rodri guez et al.ll2006h 
rim that surrounds chondrules and the interstitial matrix ma- 
terial. Figure [12] however, does not show this fine structure 
as the physical processes responsible for it were not modeled. 
Yet, the observed fact that chondrules in chondrites are sep- 
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arated by dust and not clustered together supports the main 
idea of this paper: that dust is accreted by chondrule-sized 
and perhaps somewhat larger particles (compounds) but not, 
e.g., by planetesimals. How then did the rim-matrix distinc- 
tion originate? Two scenarios can be envisioned: i) a period 
of nebula dust sweep-up and compaction; or ii) shock waves 
in the parent body. 

The first scenario concerns a moderately intense turbulent 
environment (i.e., high a) in which chondrules are largely un- 
able to stick, so that most of the dust is accreted by individ- 
ual chondrules. These bouncing chondrules quickly compact 
each other's rims, while grazing collisions may also result in 
dust being partially stripped away or eroded from the rims. 
Presumably, a steady-state between rim accretion and erosion 
is established, where some of the dust is firmly attached to 
each chondrule and compacted, while another, more fluffy, 
component is continuously eroded off and reaccreted to the 
chondrule surfaces. Any of this latter, loosely bound phase 
which remains attached to chondrules at the point they are ac- 
creted to their parent planetesimal would be easily stripped 
away in the abrasive environment of the accreting planetesi- 
mal to become 'matrix.' 

Alternatively, a much more gentle collisional environ- 
ment may be considered in which big compounds form very 
quickly, and then continue to grow to planetesimal sizes. The 
dust is then primarily accreted by compounds, though, as we 
have argued in 34.11 a significant proportion of the dust is 
always accreted by individual chondrules. As collisional en- 
ergies stay low no identifiable rims are formed; and any rim 
signature might anyway easily abrade off on the parent plan- 
etesimal. In this scenario, the fine-grained dust rims might 
result from later processes on the parent body. Specifically, it 
has been suggested that shock waves through these planetes- 
imals (caused, for example, by viole nt collisions with other 
plane tesimal s) will compact the du s t (|T rigo-Rodrigu ez et alj 
2006). In the Trigo-Rodriguez et al.l (120061) model the highest 
compaction of the dust takes place near the solid chondrule 
surface. Thus, it is only during the planetesimal stage that 
rims become distinguishable from the matrix. 

There are a number of constraints the rim formation mecha- 
nisms must satisfy. For instance, collisions must be energetic 
enough to compact the rims significantly to explain the high 
filling factors observed in chondrites. In the nebula forma- 
tion scenario, therefore, more energetic collisions are required 
than provided by the model we present here. From the argu- 
ments given at the end of 32.3.11 we can estimate the veloc- 
ities required to compact the accretion rims into nearly ran- 
dom closely packed configurations to be 10 5 ' 2 higher than the 
sticking velocity (eq. fT4l ). or v ~ 0.1 km s _1 . Clearly, other 
compaction processes are needed than can be provided by col- 
lisions in turbulence, even between compounds approaching 
St = 1 in size. One possibility, given the preference for melt- 
ing of chondrules by M ach 7 shock waves dDesch & Connolly! 
2002: lHoodet~aT]|2005l) . is that the plausibly more numerous, 
more prevalent, weaker shocks which are consequently expe- 
rienced even more routinely by particles, would provide this 
range of collisional veloci ties for ch ondrule-size particles and 
their fractal aggregates (ICieslall2006l) . Another constraint is 
the rim-matrix distinction. In the parent body shock scenario 
it must therefore be shown that this distinction una mbiguously 
results from these shocks. iTrigo-Rodriguez et alJ (120061) pro- 
vide a qualitative idea on how this mechanism operates, and it 
would be desirable for this hypothesis to be backed-up more 
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Fig. 13. — Model predictions for the thickness of the chondrule dust rim. 
The amount of dust accreted before the chondrule's incorporation into a com- 
pound determines the dust rim mass (m r j m ). (A) Scatter plot of chondrule 
mass m c h against m r j m for 200 chondrules, chosen randomly at the conclu- 
sion of the run. Two models are shown: a = 10~ 4 and a = 10~ 6 with the 
latter population being shifted by a factor of 10 for clarity. The least-squares 
power-law fits are given by the dashed lines, which have exponents of 0.93 
and 0.91, respectively. (B) Dust-chondrule mass ratio (y-axis) at the time of 
its accretion into a compound (x-axis). 

quantitatively by, e.g., sophisticated numerical simulation. 

Yet a third observational constraint is the lin ear correlation 
between chondrule radius and rim thickness (Metzl er" et alj 
ll992t|Paque & Cuzzil 19971) . Figure Q~3]iilustrates this point in 
the context of our accretion model. In Fig. [13K the chondrule- 
mass (m c h, x-axis) is plotted against the rim mass (m^m, y- 
axis) for the two models discussed in detail in 34.11 crosses 
for the a = 10 -6 model and circles for the a = 10 4 model. 
(In Fig.[T3K the a = 10~ 6 points are arbitrarily offset verti- 
cally by a factor of 10 3 for reasons of clarity.) m^m is defined 
as dust that is accreted by individual chondrules, before they 
become incorporated into a compound. The w c h-f7J r i m relation 
is shown for 200 chondrules, randomly selected from the ini- 
tial distribution. The dashed lines show the best fit having 
slopes of 0.93 and 0.91, respectively. The near- linear trend of 
rrinm with chondrule mass i s obvious but the sp read is large, 
as seen in actual chondrites (Metzl er" et alj|1992l) . Figure [T3B 
shows the accretion history of these compounds: the mass 
ratio, mrim/ffJch, is plotted (i/-axis) against the time at which 
the chondrule is swept up by a compound. The 'a = 10 -4 
chondrules' lie to the left of the 'a = 10~ 6 chondrules,' re- 
flecting their shorter collision times. The initially linear trend 
breaks down at later times as the density of free-floating dust 
decreases. Although many processes contribute to the spread 
in the data points of Fig.|6}\ - for example, differences in ve- 
locity field (linear/square-root regime) during the simulation 
and the bouncing history of chondrules - the stochasticity in 
the chondrule-compound accretion time is the main contribu- 
tor. Note also the pile-up of particles in the a = 10 4 model 
near t ~ 10 3 yr, the final time of the simulation: these are the 
chondrules that remained single during the entire simulation. 

The relation of rim thickness with chondrule size can be 
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natur ally understood as the outcome o f a nebula accretion pro- 
cess (iMorfill et al.lll998tlCuzzHl2004l Fig.O- The observed 
linear relationship in chondrites therefore suggests this rela- 
tionship should somehow have survived further processing. 
As dust rim accretion in the violent collisional environment 
differs from the non-fragmentation environment in which our 
simulations are performed, it still remains to be shown that the 
linear relationship is maintained after fragmentation/erosion 
sets in. Alternatively, if the imprints of nebula dust-accretion 
are destroyed during parent body accretion, a different mech- 
anism must explain the observed relationship. 

Future work - e.g., experimental work on rim-chondrule 
size ratios and more advanced theoretical models - must de- 
termine which of the two scenarios described above is more 
likely. Dust fragmentation and additional compaction mech- 
anisms may be included into the present model. Increas- 
ingly energetic collisions (when compounds grow towards the 
St = 1 barrier) may disrupt compound objects without strip- 
ping the rims entirely off of individual chondrules, and in 
doing so may compact the surviving fine grained rims fur- 
ther than the c d = 0.33 limit we have adopted in this study. 
Also, size distributions in the fine-grained component might 
also allow a greate r degre e of packing than in our models and 
iBfum & SchrapleJ ([2004 ) expect, in which the grains are all 
monodisperse. 

6. SUMMARY 

We have investigated a chondrule-dust aggregation mecha- 
nism in which the fine-grained dust acts as the glue that allows 
chondrules to stick. We argue that the energy in collisions is 
sufficient to compress directly accreted material, which ini- 
tially has a porous 'fairy-castle' structure, into a more com- 
pact state having a porosity th at is roughly 67% (based on 
compaction measurements by iBlum & Schraplerj 120041 and 
theoretical arguments). We have applied this model to a vari- 
ety of questions regarding the meteorific record: the relation 
of individual chondrules to their fine-grained dust rims, the 
internal structure of the chondrites, and the ability of growth 
by sticking to surpass the meter-size barrier. This study only 
starts to address these questions; more sophisticated models 
are needed to answer detailed questions on the structure of 
the meteorites. 

We find that porous accretion rims do indeed cushion col- 
lisions and facilitates growth to compound objects contain- 
ing many rimmed chondrules, but this growth is limited to 
30 - 100 cm radius objects under the most favorable condi- 
tions. This is because the chondrule component sweeps up 
all the local dust in a short time (10 2 - 10 4 yr, depending 
on nebular location) and these compounds experience higher 
relative velocities during their growth stage. Subsequent col- 
lisions merely pack the existing rims down further, so that the 
system ultimately reaches a dead-end steady state where col- 
lisions only result in bouncing, or possibly disruption. Other 
conclusions from this study are: 

• Compound growth works best in a quiescent environ- 
ment (high gas density, low a values). In a more violent 
collisional environment (a ~ 10~ 3 " -2 ) it is difficult to 
accrete dust fractally on chondrules surfaces and the en- 
ergetic collisions between compounds quickly compact 
the remainder such that collisional growth is quickly 
terminated. 



• The importance of the other parameters on the accre- 
tion process is mostly minor. The radial location does 
not affect the final growth of the compounds, although 
timescales are longer at larger R. Ice, rather than sil- 
icate, as the sticking agent will lead to bigger com- 
pounds (but we note icy grains do not dominate the me- 
teoritic record). 

• In no single model do compounds grow to planetes- 
imal sizes. Either turbulent or systematic velocities 
are too high for the porous dissipation mechanism, the 
St = 1 size being the bottleneck; but we may have over- 
estimated the systematic velocities in this study by not 
accounting for particle collective effects in low turbu- 
lence nebulae. 

• We anticipate that the dust accreted by individual chon- 
drules - before chondrules coagulate into compounds - 
will finally end up as the chondrule rim. In strong tur- 
bulent models this fraction is very high, but it remains 
significant (tens of percents) even in the collisionally 
gentle models. 

• However, at the current state of the art of this model, 
fine-grained accretion rims have a porosity significantly 
larger than seen in actual rims. Other compaction pro- 
cesses are not hard to envision, such as higher velocity 
collisions by larger mass compound objects, or nebula 
shock waves (peripheral to those energetic enough to 
melt chondrules). These remain to be modeled. 

• When we define the rim as the dust swept up by indi- 
vidual chondrules, we find very good agreement with 
the nearly linear (average) correlation between rim 
thickness and under lying chondrule radius seen in CM 
and C V chondrites dMetzler et all 1 9921: IPaque & Cuzzil 
U9971) . 

Future work will focus on two aspects of our coagulation 
model: 

• An improvement of the collisional physics, i.e., includ- 
ing fragmentation as a collisional outcome for veloci- 
ties above ~ m s ; and a refinement of the characteri- 
zation of the compound structure, e.g., to allow dust to 
compact to higher filling factors. 

• Inclusion of a proper description of the vertical struc- 
ture of the nebula, i.e., taking account of phenomena 
such as settling, collective effects, and shear turbu- 
lence. Especially, the transition from global- (a) to 
shear-turbulence is important, and in future work we 
will address this critical junction. 
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